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PREFACE 


This  Technical  Note  (TN)  describes  the  automated  moisture  and  temperature 
forecast  model  that  currently  exists  at  the  Air  Force  Global  Weather  Central 
(AFGWC).  The  AFGWC  uses  the  forecasts  from  this  model  to  provide 
environmental  support  to  the  USAF,  the  U.S.  Army,  and  other  DOD  agencies. 

This  technical  note  provides  detailed  information  on  the  scientific  and 
computational  methods  used  by  the  model  to  satisfy  military  requirements. 

This  informaton  will  be  of  particular  interest  to  personnel  in  military  and 
non-military  units  associated  with  DOD  activities  who  use  AFGWC  cloud  and 
temperature  forecasts  in  operational  or  developmental  work. 

We  are  indebted  to  the  many  people  who  have  developed  and  maintained  the 
"cloud  models"  over  the  past  ten  years.  The  many  notes,  letters,  studies,  and 
reports  which  they  prepared  made  this  TN  possible.  Special  thanks  go  to  Major 
Arnold  L.  Friend  and  Dr.  Kenneth  B.  Mitchell. 
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1.  INTRODUCTION 


1.1  General  Motivation 


The  need  for  cloud  forecasting  models  Is  based  upon  requirements  of  the 
Department  of  Defense.  Clouds,  because  of  their  restriction  to  visibility, 
severely  hamper  Air  Force  and  Array  decision-makers  who  must  operate  in  the  air 
and  ground  environment.  The  cloud  models  at  AFGWC  exist  in  order  to  provide 
input  to  these  decision-makers.  They  must  know  whether  or  not  a  target  can 
be  acquired  visually  or  through  cloud  or  moisture-sensitive  guidance  systems. 
The  types  of  targets  range  from  the  end  of  a  runway  for  landing  aircraft  to 
battlefield  conditions  for  weapon  delivery.  Because  of  these  requirements, 
the  need  for  cloud  forecasting  exists.  The  AFGWC  cloud  models,  described  in 
the  following  sections,  are  designed  to  support  that  need. 

1.2  Philosophy 

The  Five-Layer  (5LAYER)  model  is  an  automated,  synoptic  scale,  cloud  fore 
casting  program.  Its  primary  purpose  is  to  produce  cloud  forecasts  over  two 
hemispheric  domains.  Figures  1  and  2  show  the  forecast  domain  in  the  Northern 
Hemisphere  (NH)  and  Southern  Hemisphere  (SH)  respectively.  Temperatures  and 
moisture  are  forecast  and  combined  to  produce  other  forecast  parameters  such 
as  icing  and  precipitation  type.  In  this  or  any  model,  the  user  needs  to 
understand  the  characteristics  or  properties  of  the  model.  With  this 
understanding,  the  user  will  be  able  to  anticipate  model  strengths  and 
weaknesses  in  order  to  optimize  the  application  to  operational  problems. 

1.3  History 

Automated  cloud  forecasting  at  AFGWC  began  in  the  late  19603  with  the 
introduction  of  a  trajectory  model  having  a  grid  resolution  of  approximately 
200  nm  over  a  hemispheric  domain. 

Later,  the  basic  cloud  and  temperature  forecasts  were  combined  into  other 
forecast  elements  such  as  icing  and  precipitation  amount.  In  the  early  1970s, 
the  forerunner  of  the  current  5LAYER  model  was  developed  with  a  relocatable 
grid  of  100  nm  resolution. 

Also  in  the  early  1970s,  a  high  resolution,  limited  area,  short  range 
forecast  capability  was  appended  to  the  basic  model  to  take  advantage  of  the 
Three-Dimensional  Nephanalysls  Model  ( 3DNKPH)  (Fye,  1978).  The  3DNKPH  was 
succeeded  in  August  1983  by  the  Real  Time  Nephanalysls  (RTNEPH)  Model.  An 
Improved  mathematical  technique  for  calculating  the  trajectories  was  adopted 
during  the  mid  1970s.  About  the  same  time,  the  100  nm  grid  resolution  was 
applied  over  the  entire  hemispheric  domain  thereby  eliminating  the  need  for  a 
relocatable  grid. 
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Fig.  1.  The  Whole-mesh  and  Half- mesh  Octagons  for  the  Northern 
Hemisphere.  Tha  domain  of  the  Reference  Grids  (solid  outer 
border)  Is  also  plotted.  The  indices  (Iy,Jy)  and  (I2,J2) 
designate  the  coordinates  for  the  whole-mesh  and  half  mesh 
versions,  respectively,  with  respect  to  the  Octagons  and  the 
Reference  Grids.  The  dashed  lines  indicate  that  data  on  the 
Octagons  are  actually  stored  In  a  rectangular  database.  The 
whole  mesh  grid  spacing  is  displayed  in  the  upper  left  corner. 
From  Hoke  et  al .  (1985). 
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Fig.  2.  Same  as  Fig.  1  except  for  the  Southern  Hemisphere 


In  1986,  a  new  initialization  module  was  Implemented.  Thla  allowed  the 
FLAYER  to  be  Initialized  directly  from  Lhe  RTNEPH  database.  Another  major 
change  was  the  ability  of  5LAYER  to  use  Its  forecast  to  Initialize  the 
moisture  field  If  the  a-.alysls  data  was  older  than  a  specified  threshold. 

This  threshold  Is  adjustable  and  presently  Is  set  at  4  hours. 

1.4  Current  Specifications 

1.4.1  Forecast  Elements 

The  SLAYER  model  forecasts  cloud  and  temperature  In  three-hour  increments 
out  to  48  hours  in  the  NH  and  24  hours  in  the  SH.  The  choice  of  the 
three-hour  time  step  is  based  on  convenience  rather  than  necessity,  Wind 
components  obtained  from  the  Air  Weather  Service  Global  Spectral  Model  (AWS 
GSM)(Stobie,  1986),  40  wave  and  12  level,  are  available  In  three-hour 
Increments;  therefore,  trajectories  are  computed  to  represent  parcel  movement 
over  that  same  time  Increment.  In  this  approach,  the  length  of  the  time  step 
Is  unrelated  to  any  numerical  stability  criteria  because  the  Lagranglan 
advectlon  scheme  of  5LAYKR  does  not  explicitly  solve  finite  difference 
approximations  to  differential  equations. 

Temperature  forecasts  In  combination  with  cloud  forecasts  can  produce 
additional  meteorological  elements.  Air  Force  and  Army  decision-makers  are 
not  only  concerned  with  clouds  but  also  with  any  meteorological  element  that 
effects  the  accomplishment  of  their  mission.  Two  Important  elements  are 
precipitation  and  Icing.  Conditions  that  produce  aircraft  Icing  must  be 
anticipated  by  the  operational  planner.  The  5LAYER  model  uses  several  methods 
to  produce  forecasts  of  precipitation  amount,  precipitation  type,  and  icing. 
Dew  point  depressions,  static  stability,  and  cloud  types  are  forecast  to 
determine  these  elements. 

The  discontinuous  nature  of  cloud  coverage  suggests  the  need  to  forecast 
clouds  on  a  small  resolution  grid  system.  However,  constraints  such  as 
computer  capabilities  and  the  need  to  meet  operational  time  requirements  cause 
the  spatial  resolution  to  be  less  than  ideal.  The  5LAYBR  model  forecasts  for 
five  layers  in  the  vertical  and  on  grid  points  separated  by  100  run.  Even 
though  the  horizontal  grid  spacing  is  approximately  equal  to  the  smallest  grid 
spacing  of  many  operational  dynamic  models,  this  spacing  causes  considerable 
dilution  of  the  available  Information.  Meteorological  satellites  routinely 
resolve  cloud  elements  to  one- third  of  a  mile.  Averaging  or  compacting  this 
Information  to  a  grid  resolution  of  100  nm  causes  smoothing  of  the  data.  The 
problem  of  data  compaction  will  be  discussed  later.  Users  of  cloud  model 
forecasts  should  constantly  remind  themselves  that  cloud  forecasts  represent 
average  cloud  conditions  over  a  large  volume.  This  same  volume  is  the 
"parcel”  of  air  that  is  referred  to  throughout  this  memo. 
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1.4.2  Forecast  Mechanism 

Generally,  SLAYER  uses  a  quasi-  Lagranglan  advectlon  scheme  to  determine 
the  advected  elements.  A  pure  Lagranglan  scheme  follows  a  parcel  throughout 
the  forecast  period.  SLAYER  follows  the  parcel  In  three  hour  increments  until 
the  end  of  the  forecast  period.  A  characteristic  of  this  scheme  Is  that  the 
trajectory  remains  fixed  at  one  end  while  the  other  end  specifies  the  source 
of  the  parcel.  This  upstream  trajectory  requires  a  grldded  analysis  to 
provide  values  of  the  element  at  the  source  point.  The  value  of  the  element 
Is  determined  by  Interpolation  to  the  source  point  and  modified  as  it 
traverses  the  trajectory  path.  Since  all  trajectories  terminate  at  a  grid 
point,  a  grldded  forecast  field  results.  This  procedure  Is  repeated  with  each 
forecast  field  serving  as  the  Initial  field  for  subsequent  forecast  Increments 
until  the  desired  forecast  length  Is  reached. 

The  motivation  to  use  this  type  of  forecast  scheme  lies  In  the 
characteristically  discontinuous  nature  of  clouds.  This  discontinuous  nature 
Is  not  easily  adaptable  to  conventional  advectlon  schemes  which  require  a 
smooth,  continuous  variation  of  the  advected  element.  Finite  difference 
approximations  to  differential  equations  perform  best  when  gradients  of  the 
element  can  be  accurately  determined.  The  quasi- Lagranglan  approach  Is  an 
acceptable  alternative  since  gradients  of  clouds  are  not  resolvable  by  current 
numerical  models. 

1.5  Overview 

The  remaining  sections  of  this  technical  note  describe  specific  features 
of  the  cloud  forecasting  models  at  APGWC.  Section  2  details  the  Initialization 
procedures  of  the  model  for  both  moisture  and  temperature.  Section  3 
discusses  the  computation  of  trajectories  used  to  advect  the  moisture  and 
temperature.  Section  4  describes  how  other  physical  processes  are  combined 
with  the  advectlon  calculation  to  obtain  the  moisture  and  temperature 
forecasts.  The  section  then  goes  on  to  explain  how  several  additional 
quantities  of  meteorological  interest  are  derived  from  the  moisture  and 
temperature  forecasts.  Section  5  describes  the  procedures  used  to  obtain  a 
limited  area  cloud  forecast  of  up  to  9  hours  at  25  run  resolution.  Section  6 
discusses  the  procedures  used  to  forecast  clouds  in  the  tropics.  In  Section  7 
verification  of  the  model  is  reviewed.  Several  appendices  are  Included  to 
give  explicit  details  of  some  of  the  more  complex  mathematical  formulations 
and  the  data  used  for  the  moisture  to  cloud  amount  conversions. 
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INITIALIZATION 


2.1  General  Characteristics 

To  produce  a  numerical  lorecast,  one  must  first  construct  analysis  fields 
of  the  meteorological  elements  to  be  forecast.  Two  such  elements  of  major 
interest  at  AFGWC  are  temperature  and  moisture.  Because  of  its  world-w' le 
mission,  AFGWC  must  forecast  these  elements  on  a  global  basis.  To  supp  ?t 
this  requirement,  automated,  hemispheric,  upper  air  analysis  models  produce 
temperature  and  dew-point  depression  analyses.  In  addition,  an  automated  cloud 
analysis  model  supplies  initial  cloud  fields  by  merging  visual  and  Infrared 
(IR)  satellite  data  with  conventional  meteorological  data.  The  relationship 
between  analysis  and  forecast  fields  requires  detailed  discussion  of  cloud  and 
temperature  Initialization  to  gain  a  better  understanding  of  not  only  the 
initialization  procedures,  but  also  the  quality  and  verification  of  forecast 
variables.  A  conceptual  flow  chart.  Figure  3,  summarizes  the  Initialization 
procedures  that  will  be  discussed  In  the  following  sections. 

2.1.1  Motivation 

The  quality  of  a  forecast  is  no  greater  than  the  quality  of  the  analysis 
upon  which  it  is  based.  For  thi3  reason,  the  Initialization  of  the  cloud 
forecast  models  is  designed  to  be  as  consistent  and  accurate  as  possible.  The 
input  for  the  cloud  forecast  models  comet  from  the  RTNEPH.  The  automated 
RTNEPH  analyses  contain  Information  on  the  amount  of  cloud  present.  Its  base 
and  top  for  up  to  four  cloud  layers.  This  Information  Is  based  on  satellite 
analyses  and  conventional  meteorological  data  (surface  observations, 
radiosondes,  and  aircraft  reports).  Many  factors  limit  the  consistency  and 
accuracy  of  the  initialization  of  the  moisture  and  temperature  variables. 

These  Include  the  grid  resolution  of  Ihe  model,  the  timeliness  of  the 
observation,  and  our  understanding  of  the  physics  of  the  atmosphere  as  related 
to  clouds  and  moisture. 

The  grid  resolution  affects  the  initialization  of  moisture  because  of  data 
compaction.  "Compaction”  refers  to  the  process  of  averaging  data  at  one 
resolution  to  yield  data  at  a  coarser  resolution.  Very  high  resolution 
satellite  data  Is  compacted  to  25  nm  resolution  by  the  RTNEPH  and  then  further 
compacted  to  100  run  to  initialize  the  forecast  model.  Each  reduction  In 
resolution  tends  to  decrease  the  accuracy  of  the  analysis  because  the 
averaging  process  decreases  the  sharp  distinction  between  moist  and  dry 
areas.  Consistency  between  the  analysis  of  moisture  by  the  RTNEPH  and  the 
Initialization  fields  of  the  forecast  model  is  difficult  because  the  analysis 
is  predominantly  based  on  total  cloud  data  from  meteorological  satellites 
while  the  initialization  of  the  forecast  models  must  be  on  their  distinct 
pressure  surfaces. 
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Since  the  satellite  data  used  In  the  HTNEPH  Is  from  polar  orbiting 
satellites  and  the  5LAYER  grid  covers  the  extratroplcal  regions  of  the 
Northern  and  Southern  hemispheres,  the  data  used  to  analyze  the  cloud  field 
can  be  8  or  more  hours  old.  Instead  of  using  data  this  old,  the  model  will 
Initialize  from  cloud  forecasts  valid  at  the  current  cycle  time.  Even  with 
forecast  data,  however,  fast-moving  cloud  systems  create  discontinuities  near 
the  edges  of  overlapping  cloud  coverage.  The  discontinuities.  In  turn,  cause 
inconsistencies  in  the  initialization. 

In  order  to  convert  from  the  cloud  elements  sensed  by  the  satellites  to 
the  moisture  used  as  the  forecast  element,  empirical  relations  have  been 
developed.  While  these  relations  have  proven  to  be  usable,  they  still  fall 
short  of  the  accuracy  which  would  be  obtained  by  a  specific  physical 
relationship. 

The  initialization  of  temperature  has  the  inverse  problem  of  grid 
resolution  since  the  analyzed  temperatures  are  based  upon  whole-mesh  (200  nm) 
temperatures  which  are  derived  from  the  High  Resolution  Analysis  System 
( HI HAS )  which  is  based  on  a  resolution  of  2.5  degrees  of  latitude.  The 
initial  temperature  values  are  obtained  by  interpolating  to  the  100  nm 
resolution  of  the  cloud  model. 

2.1.2  Requirements 

The  sole  purpose  of  the  initialization  Is  to  provide  an  accurate  and 
consistent  set  of  moisture  and  temperature  values  with  which  to  start  the 
forecast  at  each  model  grid  point.  As  such,  the  moisture  values  on  the 
various  model  surfaces  must  be  consistent  enough  to  be  recombined  into  a  total 
amount  that  Is  representative  of  the  total  cloud  at  that  grid  resolution. 
Temperature  must  also  satisfy  certain  conditions.  Specifically,  the  lapse 
rate  must  be  constrained  to  avoid  superadiabatlc  conditions.  Additionally, 
model  surfaces  which  Intersect  the  terrain,  which  can  often  be  at  850  mbs, 
must  have  the  below  ground  grid  points  identified. 

<1.2  Cloud  Initialization 

At  AFGWC ,  initial  cloud  fields  are  constructed  using  the  Real-Time 
Nephanalysis  Model  (HTNEPH).  This  model  combines  satellite  data,  rawlnsonde 
data,  conventional  surface  observations,  and  aircraft  reports  to  provide  up  to 
four  layers  of  cloud  information  (in  percent  cloudiness)  on  a  grid  spacing  of 
approximately  25  nm.  However,  because  of  several  operational  considerations, 
this  highly  detailed  information  must  be  compacted  to  a  coarser  grid  spacing. 
The  AFGWC  grid  system  for  cloud  and  lemperature  Initialization  has  a 
horizontal  spacing  of  approximately  100  nm  and  a  vertical  resolution  of  five 
layers . 
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2.2.1  Data  Compaction. 

The  AFGWC  200  nm  grid,  subsequently  Identified  as  whole  mesh,  is  a  subset 
of  the  RTN8PH  horizontal  grid.  The  major  difference  between  these  two  grid 
systems  is  that  the  RTNEPH  grid  has  points  every  25  nm  and  is  referred  to  as 
eighth  mesh.  The  more  detailed  resolution  consists  of  262,144  (4096  x  64) 

horizontal  points  per  hemisphere.  At  each  of  these  points  the  RTNEPH  analyzes 

up  to  four  layers  of  cloud  information.  At  eighth- mesh  on  a  global  scale,  the 
5LAYBR  model  would  need  to  initialize  2,621,440  (262,144  grid  points  per 
hemisphere  x  2  hemispheres  x  5  levels)  grid  points.  The  currently  available 
computer  hardware,  coupled  with  the  time  constraints  at  AFGWC,  make  it 
impossible  for  a  prediction  model  to  treat  each  of  these  points. 

The  vertical  structure  of  the  RTNEPH  grid  consists  of  four  floating  cloud 
layers  at  each  point.  At  any  particular  point,  some  or  all  of  the  layers  may 
be  cloud  free.  For  each  layer,  the  cloud  top,  bottom,  amount,  and  type  are 
given.  The  layers  are  sorted  based  on  cloud  tops,  from  highest  to  lowest. 
Total  cloud  (in  percent)  is  also  given  for  each  point,  as  is  a  time  flag.  For 

5LAYER  to  use  the  RTNEPH  data,  it  must  be  spread  to  5  vertical  layers  and 

compacted  to  a  half-mesh  horizontal  grid. 

The  RTNEPH  data  are  compacted  to  reduce  the  volume  of  data,  while 
retaining  most  of  the  larger-scale  information.  For  example,  to  construct  a 
value  at  a  half-mesh  (100  nm)  grid  point,  the  data  of  the  25  RTNEPH  grid 
points  (in  the  eighth-mesh  grid  weighting  system)  that  surround  this  point  is 
used.  This  data  compaction  of  eighth-mesh  points  is  weighted  so  that  each 
eighth-mesh  grid  point  influences  the  final  compacted  product  equally.  The 
scheme  in  Figure  4  is  applicable  for  all  half-mesh  grid  points  that  are 
internal  to  the  grid  borders.  For  those  half-mesh  points  coincident  with  the 
grid  border  the  weighting  scheme  is  changed  so  that  the  number  of  eighth  me3h 
points  considered  in  the  half-mesh  is  correspondingly  reduced. 

To  spread  the  clouds  vertically,  the  top  and  base  of  each  layer  is  checked 
in  turn  to  see  which  5 LAYER  fixed  layers  they  overlap.  The  RTNEPH  layer 
amount  is  inserted  into  the  overlapped  5 LAYER  layers  (unless  the  receiving 
layer  already  has  a  larger  amount  in  it).  The  process  is  repeated  for  each 
RTNEPH  layer  that  has  cloud  in  it.  Empty  layers  are  skipped. 

The  vertically  spread  data  is  then  multiplied  by  a  weighting  factor  and 
added  to  the  adjacent  half -mesh  points.  It  can  be  added  to  up  to  four  points, 
depending  on  its  position  relative  to  those  points.  The  weighting  factor  Ls 
also  added  to  a  weight  holding  array.  The  factor  is  chosen  so  that  no  matter 
how  many  half- mesh  points  an  eighth-mesh  point  adds  into,  its  overall 
contribution  is  the  same.  This  data  composition  of  eighth-mesh  points  is 
weighted  so  that  each  eighth- mesh  grid  point  influences  the  final  compacted 
product  equally. 

After  the  weighted  eighth- mesh  data  is  added  to  the  half-mesh  octagon,  the 
half-mesh  values  are  divided  by  the  respective  sum  of  weights.  The  resulting 
point  layer  cloud  amounts  are  adjusted  so  the  point  total  cloud  ls  within  ?% 
of  the  weighted  average  RTNEPH  total  cloud  at  the  point. 


HORIZONTAL  COMPACTION 


WEIGHTS  ASSIGNED 


INTERIOR 

1.00 

SIDE 

0.50 

OOFNER 

0.25 

Fig.  4.  Weights  assigned  to  elghth-mesh  (x)  grid  point  values 
while  compacting  to  half-mesh  (H).  The  locations  (Interior, 
side,  corner)  are  relative  to  the  flve-by-flve  box  of  elghth-mesh 
grid  points. 
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2.2.2  Total  Cloud  Computation 

The  layered  cloud  Initialization  described  above  Is  necessary  to  produce 
forecasts  of  layered  clouds.  Although  total  cloud  Is  not  a  direct  forecast 
parameter,  It  Is  inferred  for  verification  purposes  and  some  applications. 

The  verification  relationship  between  cloud  forecasts  and  cloud  analysis  must 
be  made  via  total  cloud  because  observations  of  layered  cloud  are  not  frequent 
or  detailed  enough  to  satisfy  verification  requirements.  Thus,  a  precise 
methodology  Is  required  to  establish  the  relationship  between  layered  and 
total  cloud.  This  methodology  must  produce  a  total  cloud  on  the  half  mesh 
grid  which  closely  matches  the  total  cloud  amount  obtained  by  horizontal 
compaction  of  the  RTNEPH  total  cloud.  The  total  cloud  values  of  the  H'FNEPIl 
Influence  the  layered  cloud  Initial  values  of  the  forecast  model  and  provide 
verification  of  previous  forecasts. 

The  total  cloud  at  a  point  is  computed  from  the  statistical  union  of  layer 
cloud  amounts,  maximum  layer  amount  and  average  separation  of  the  cloud 
layers.  Since  layer  cloud  may  be  separated  by  thousands  of  meters,  an 
estimate  of  how  well  these  layers  are  correlated  Is  needed. 

For  example,  assume  each  of  five  cloud  layers  contain  50  percent  cloud. 

If  the  layers  are  perfectly  correlated  (r  -  1),  the  total  cloud  amount  equals 
50  percent,  the  maximum  layer  amount.  If  the  layers  are  independent  (r  -  0), 
the  total  cloud  amount  is  about  97  percent,  the  statistical  union  of  the 
Individual  layers. 

The  actual  layer  correlation  is  somewhere  between  zero  and  one,  so  the 
total  cloud  In  the  above  example  would  really  be  between  50  and  97  percent. 

The  problem  Is  how  to  estimate  this  correlation.  A  conceptual  illustration  of 
the  vertical  cloud  stacking  problem  Is  given  In  Figure  5. 

A  factor,  based  on  layer  separation,  is  used  to  estimate  the  layer 
correlation.  The  basic  assumption  is  that  layers  separated  by  the  depth  of 
the  troposphere  (assumed  to  be  11000  meters)  are  uncorreiated  (r  *  0),  while 
those  layers  immediately  adjacent  to  each  other  are  perfectly  correlated  (r 
1).  Intermediately  separated  layers  have  correlations  which  are  assumed  to 
vary  linearly  with  mean  layer  separation: 

Correlation  =  1  -  (Helghtlayer  1  -  Helghtlayer  2  )/  11000  (1) 

For  more  than  two  cloud  layers,  the  average  of  all  possible  pair-wise 
separations  of  the  cloud  layers  is  divided  by  11000  and  then  subtracted  from 
one  to  estimate  the  correlation. 

The  final  algorithm  used  to  obtain  the  total  cloud  is: 

(1)  Determine  the  maximum  layer  cloud  amount.  This  Is  equal  to  the 
minimum  total  cloud. 
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Side  and  top  view  of  vertical  cloud  stacking  problem 


(II)  Calculate  the  statistical  union  of  the  layer  cloud  amounts.  This  Is 
equal  to  the  maximum  total  cloud. 

(III)  Calculate  the  mean  separation  of  the  cloudy  layers,  divide  by 
11000,  and  subtract  from  1.  This  estimates  the  correlation. 

NOTE:  If  there  are  n  layers,  there  are  n«(n-l)/2  possible  unique 
pairwise  layer  separations  which  need  to  be  averaged  to  obtain  the  mean 
separation. 

(lv)  Use  results  from  (1)  -  (111)  to  estimate  the  total  cloud: 

Total  =  Min  Total  ♦  (Max  Total  -  Min  Total)  "  Correlation  (2) 

This  computed  total  Is  then  used  to  adjust  the  layers  (Figure  6)  to  make 
them  match  reality  (assumed  to  be  the  RTNBPH  derived  cloud  total).  The 
adjustment  factor  Is  the  RTNBPH  total  divided  by  the  computed  total.  If  the 
adjustment  factor  Is  greater  than  1.02  or  less  than  .98,  all  the  layers  are 
multiplied  by  the  adjustment  factor  to  Increase  or  decrease  the  layer  totals. 
Then  the  computed  total  cloud  Is  recalculated  (as  above)  and  the  adjustment 
repeated.  The  layers  will  be  adjusted  up  to  three  times,  or  until  the 
adjustment  factor  is  In  the  range  .98  -  1.02.  In  either  case,  the  total  cloud 
used  to  initialize  the  5LAYBR  database  Is  the  RTNBPH  total  cloud. 

2.2.3  Moisture  Determination. 

2. 2.3.1.  Moisture  In  Cloudy  Areas. 

Condensation  pressure  spread  (CPS)  Is  the  moisture  parameter  used  by 
5LAYER  to  forecast  cloud  amount.  CPS  is  the  difference  (In  millibars)  between 
the  pressure  of  an  air  parcel  and  the  pressure  at  which  condensation  takes 
place  If  the  parcel  Is  lifted  dry  adlabatically .  Written  mathematically, 

CPS  =  P  -  P„  (3) 

where  P  Is  the  pressure  of  the  parcel  and  P8  Is  the  saturation  pressure. 

See  Appendix  A  for  the  derivation  of  equation  (3).  To  calculate  CPS  at  a 
fixed  pressure  level,  P,  only  Pa  must  be  specified.  An  approximate, 
simplified  relationship  between  CPS  In  mb  and  dew-point  depression  (DPD)  In  K, 
can  be  established  as: 
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Kig.  6.  Logic  flow  for  adjusting  the  layer  cloud  amount  spread 
to  the  five  layers  in  5LAY8R  to  make  the  computed  total  cloud 
amount  near  the  observed  total  cloud  amount. 


CPS  -  [  B0(P)  +  B^P)  *  DPD  )  *  DPD 


(4) 


where 

Bq(P)  -  BOO  ♦  B01  *  p 
Bi(P)  B10  +  Bn  *  P 


and. 

Bqq  =  1.41985 
B01  =  1.34466E  2 

B10  =  -1.39131E-2 
Bn  =  -  6.69419E  5 

The  resulting  error  Is  less  than  4%,  when  compared  with  an  exact  conversion. 

Now  that  the  relationship  between  CPS  and  dew-point  depression  has  been 
established,  It  remains  to  relate  CPS  to  cloud  amount.  Intuitively,  small 
values  of  CPS  (small  dew-point  depressions)  correspond  to  large  values  of 
cloud  amount,  and,  conversely,  large  CPS  values  (large  dew  point  depressions) 
correspond  to  small  cloud  amounts.  Edson  (1965)  presented  empirical  curves 
relating  these  two  variables.  These  empirical  curves  are  shown  In  Figure  7 
and  tabulated  In  Appendix  B. 

There  are  several  advantages  to  using  CPS  as  the  moisture  parameter. 

First,  as  previously  discussed,  CPS  provides  the  link  between  dew-point 
depressions  and  cloud  amounts  required  for  Initialization  In  cloud  and 
cloud-free  areas.  Second,  CPS  also  links  changes  In  cloud  amount  to  changes 
In  vertical  motion.  For  example,  should  an  unsaturated  parcel  of  air  be 
lifted,  adiabatic  cooling  takes  place.  Simultaneously,  the  dew-point 
depression  narrows,  and  additional  clouds  form.  For  descending  motion,  the 
air  becomes  drier,  and  clouds  dissipate.  Therefore,  cloud  amounts  In  terms  of 
CPS  (units  of  pressure)  can  be  modified  directly  by  the  net  vertical 
displacement  (units  of  pressure)  of  a  parcel  of  air. 

2. 2. 3. 2  Moisture  In  clear  areas. 

Initialization  schemes  up  to  this  point  have  determined  CPS  In  areas  where 
clouds  exist.  However,  CPS  values  have  not  been  adequately  described  in 
cloud- free  regions.  For  example,  a  volume  of  air  with  no  cloud  may  be  very 
dry  or  It  may  be  so  moist  as  to  require  only  minimum  adiabatic  cooling  to 
begin  cloud  formation.  The  RTNSPH  cloud  Information  cannot  define  moisture  in 
cloud  free  regions.  Consequently,  for  these  regions  we  turn  to  an  alternate 
data  source. 
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To  provide  this  Information,  a  dew  point  depression  analysis  from  the 
AFGWC  Multi-level  Analysis  Model  (Tarbell  and  Hoke,  1979)  Is  used.  The 
Initialization  scheme  creates  the  half  mesh  (100  run)  field  by  Interpolation 
from  whole  mesh  values  as  described  In  Appendix  C. 

This  three-step  Initialization  process  Is  summarized  by  the  flow  diagram 
In  Figure  3.  First,  the  half  mesh  analysis  Is  constructed  using  RTNKPH  values 
In  cloudy  areas.  Secondly,  the  dew-point  analysis  In  cloud  free  regions  la 
determined  from  the  conventionally  analyzed  dew-point  depression  values.  Then 
a  check  Is  made  to  ensure  that  the  drier  of  the  analyzed  dew  point  depressions 
or  cloud-to-C?S  conversion  Is  used. 

However,  a  problem  has  been  created  in  that  the  moisture  field  In  the 
Initialization  scheme  consists  of  both  cloud  Information  (In  percent  cloud 
average)  and  dew-point  depressions  (In  degrees).  A  common  moisture  parameter 
needs  to  be  specified.  This  parameter,  condene**tlon  pressure  spread  (CPS), 
links  cloud  amount  and  dew-point  depression. 

2.3  Temperature  Initialization 

If  temperature,  or  any  parameter  requiring  It,  Is  forecast,  an  Initial 
temperature  field  Is  required.  At  APGWC  this  field  Is  created  by  the  High 
Resolution  Analysis  System  (HIRAS).  The  HIRAS  derives  temperatures  from  the 
geopotentlal  thickness  between  two  pressure  levels  at  six  levels  (1000  mb,  850 
mb,  700  mb,  500  mb,  A00  mb,  300  mb).  These  temperatures  are  converted  to  a 
whole  mesh  (200  nm)  grid  spacing. 

Some  modification  of  this  Information  Is  required  to  format  It  in  the  five 
level,  half-mesh  (100  nm)  grid  system  of  5L.AYRR.  Whole-mesh  temperatures  are 
horizontally  Interpolated  to  half-mesh  In  the  same  manner  as  dew  point 
depressions  (recall  section  2. 2. 3. 2.)  The  four  required  pressure  levels 
correspond  to  pressure  level  Information  from  the  analysis  model  directly. 
However,  Initialization  of  temperature  on  the  terrain-following  (gradient) 
level  requires  additional  Information  extraction  techniques.  No  attempt  Is 
made  to  use  surface  temperature  values.  Rather,  gradient  level  temperatures 
are  assigned  based  on  standard  pressure  level  Information. 

The  gradient  level  is  Initialized  by  height  Interpolation  between  fixed 
pressure  levels.  This  technique  ensures  a  representative  value  of  temperature 
Is  assigned.  Gradient  temperatures  are  thus  relatively  free  from  local, 
small  scale  temperature  effects. 
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3.  TRAJECTORY  COMPUTATION 
3.1  General  Motivation 

Cloud  and  temperature  forecast  models  at  APGWC  are  for  the  moat  part 
diagnostic  in  nature.  However,  dynamic  properties  that  do  exist  are  found 
Implicitly  in  the  trajectory  computation.  Wind  forecasts  from  the  GSM  are 
usol  to  construct  three  dimensional  trajectories  to  determine  the  path  of  an 
air  parcel.  The  input  data  for  trajectory  computation  consists  of  GSM-derived 
whole  mesh  wind  forecasts  at  three-hour  intervals  for  the  1000  mb,  850  mb,  700 
mb,  500  mb,  and  300  mb  pressure  levels. 

By  averaging  the  wind  forecasts  in  time,  the  trajectories  are  computed 
backward  from  selected  terminal  points  (grid  points)  to  origin  points.  Since 
these  trajectories  determine  the  origin  of  an  air  parcel,  they  are  computed  by 
starting  at  the  final  forecast  point.  The  model  calculates  an  "upstream" 
trajectory  in  three  hour  intervals  until  the  parcel's  origin  is  determined. 
From  these  points,  analyzed  atmospheric  elements  such  as  clouds  and 
temperatures  can  be  displaced  and  modified  as  each  parcel  traverses  its 
trajectory  path.  Figure  8  shows  a  conceptual  flow  chart  of  the  trajectory 
computation  process. 

Very  critical  to  the  accuracy  of  these  origin  points  la  the  choice  of 
approximations  to  the  wind  field.  Using  different  geostrophlc  or 
stream  function  approximations,  DJurlc  (1961)  showed  that  in  only  12  hours  the 
average  displacement  error  of  the  origin  point  ranges  up  to  20%  of  the 
Lrajectory  length.  Therefore,  the  forecast  accuracy  of  any  meteorological 
element  (clouds  and  temperatures)  is  extremely  dependent  upon  the  accuracy  of 
th<-  wind  model.  In  an  extreme  caje,  it  is  conceivable  that  the  movement  of  a 
trough  or  ridge  pattern  could  be  so  fast  or  so  slow  as  to  cause  the  forecast 
elements  to  be  180°  out  of  phase. 

Experiments  wore  conducted  in  1974  at  AFGWC  to  quantify  the  effect  of  wind 
model  errors  on  forecast  accuracy.  Because  trajectory  origin  points  are 
difficult  to  verify,  the  experiments  were  designed  to  verify  the  accuracy  of  a 
forecast  element  such  as  clouds.  One  particular  experiment  compared  cloud 
forecasts  which  were  produced  from  three  wind  models.  One  was  a  baroclinic, 
quasi  geostrophlc  model,  another  was  a  primitive  equation  (PE)  model  and  the 
thirl  was  a  "perfect  prog"  model  that  uses  successive  wind  analyses.  Overall, 
the  cloud  verification  showed  that  the  wind  model  accounts  for  25%  of  the 
standard  deviation  of  error  (STDE).  The  other  75%  was  attributable  to 
grid  associated  shortcomings,  variance  attributable  to  wind  and  cloud 
analyses,  and  cloud  forecasting  algorithms.  Of  particular  note  was  the 
superior  performance  of  the  PR  model  over  the  quasi-geostrophlc  model.  The 
STDK’s  for  the  PE  driven  cloud  forecasts  showed  a  25%  improvement.  Therefore, 
it  was  desirable  to  link  the  trajectory  computations  to  the  AWSPE  model  which 
provided  the  best  operational  wind  forecasts  at  AFGWC.  The  GSM  succeeded  the 
AWSPE  in  October  1985.  With  the  advent  of  larger  computers  and  more 
sophisticated  wind  models,  the  future  offers  great  potential  to  increase  not 
only  trajectory  accuracy  but  also  forecast  element  accuracy. 
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Kip;.  8.  Flow  of  3tep8  taken  to  compute  5LAYBR  trajectories  from 
forecasts  of  GSM  wind  components. 


3.2  Pressure- following  Levels 


To  compute  trajectories,  the  model  uses  an  unpublished  technique  developed 
by  Mr  Ralph  Jones  (currently  at  the  National  Meteorological  Center).  The 
method  Involves  using  a  Taylor  series  expansion  about  a  grid  point  to  compute 
representative  wind  components  u,  v,  andCJ.  For  the  u  component,  the  equation 
takes  the  form 


3u  9  u  9u 

u  =  u  4-  —  —  dx  +  — —  dy  +  ■ —  dP  (5) 

s  e  3x  3y  3  3P 

where  ue  and  ua  refer  to  the  u  components  at  the  grid  point  and  trajectory 
origin,  respectively.  Derivatives  on  the  right-hand  side  are  approximately 
three  hour  average  changes.  Dy  assuming 

u  -  (ue  +  ua)/2  (6) 

and  substituting  (f>)  into  (6)  for  ua,  one  obtains 

_  .  9u  3u  3u 

u  =  u  +  —  -  dx  +  -7—  dy  +  dP] 

e  2  dx  dy  dl1 

Trajectory  displacement  components  are  then  computed  from 


Ax  =  -uAt 


(8] 


The  advantage  of  using  the  partial  Taylor  series  method  is  that  trajectory 
displacements  represent  curvature  in  the  flow  more  adequately  than  iteration 
techniques.  Appendix  D  gives  a  detailed,  explicit  solution  to  (7),  as  well  as 
comparable  derivations  for  the  v  and  CJ  wind  components. 

3.3  Gradient  Level 

Trajectory  displacements  are  also  computed  in  the  terrain- following 
gradient  level.  Computations  here  are  basically  the  same  as  those  that  are 
described  in  section  3.2.  However,  the  evaluation  of  the  derivatives  in  (7) 
at  this  level  needs  further  explanation.  Also,  because  the  gradient  level  Is 
a  terrain  following  level,  special  terrain  adjustment  features  are  inherent  to 
the  trajectory  calculations  and  need  to  be  explained. 
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In  order  to  evaluate  the  horizontal  and  vertical  derivatives  In  (7)  with 
respect  to  the  gradient  level,  a  wind  field  at  this  level  is  required.  Since 
such  a  field  does  not  normally  exist,  It  Is  created  by  linear  vertical 
Interpolation  between  known  wind  components  above,  and  terrain  level  (surface) 
components  below  the  gradient  level.  The  components  above  the  gradient  level 
come  directly  from  the  next  pressure  level.  However,  the  surface- level  winds 
may  or  may  not  be  used  In  a  direct  manner.  This  decision  is  based  on  whether 
the  terrain  pressure  tieight  Is  above  or  below  1000  mb.  If  the  pressure  height 
Is  equal  to  or  greater  than  1000  mb,  the  1000  mb  wind  components  from  the 
GSM  are  used  directly.  However,  If  It  is  less  than  1000  mb,  then  the  two 
closest  pressure  levels  are  vertically  interpolated  to  calculate  u^  and 
Vfc.  In  either  case,  these  two  components  are  modified  to  account  for  surface 
friction.  The  velocities  are  modified  by  a  roughness  turning  angle  ranging 
from  eight  degrees  over  water  to  20  degrees  over  rough  terrain. 


The  vertical,  terrain  level,  wind  component,  wt ,  is  needed  for  vertical 
Interpolation  to  the  gradient  level.  Tills  terrain  induced  motion  is 
approximated  by: 


-100(u 


Dli 


9x  VL  Dy 


(9) 


where  h  Is  the  terrain  height  In  meters.  The  100  term  in  (9)  Is  an 
approximation  to  convert  the  vertical  velocity  from  10*m/s  to  lO^mb/s. 

This  scaling  preserves  consistency  with  vertical  wind  components  at  the 
pressure  levels.  Now  that  all  wind  components  above  and  below  the  gradient 
level  are  known,  ug,  vg,  and  w  g  are  determined  by  vertical 

Interpolation.  Equation  (14)  uses  these  newly  derived  components  to  produce 
terrain  following  trajectories. 


3.4  Boundary  Characteristics 
3.4.1  Terrain  Effects 


Because  of  Its  proximity  to  the  surface,  gradient  level  wind  components 
may  be  modified  to  account  for  steep  terrain  conditions.  This  terrain  field  is 
extracted  from  the  RTNKPH,  but  It  Is  modified  to  a  half  mesh  (100  nm)  format 
by  using  a  25- point  unweighted  average. 


Modification  of  the  gradient  level  trajectories  for  terrain  effects  13 
done  for  two  practical  reasons.  The  first  Is  to  prevent  an  upstream 
trajectory  from  backing  Into  the  terrain.  This  limitation  avoids  calculation 
of  trajectories  that  have  origins  which  are  below  the  terrain  surface.  At  the 
conclusion  of  each  time  step,  the  pressure  of  the  trajectory  origin  Is 
compared  to  the  Interpolated  terrain  pressure.  If  the  origin  Intersects  or  Is 
below  the  terrain  level,  the  vertical  component  of  the  trajectory  Is  adjusted 
so  that  it  Is  exactly  20  mb  above  the  terrain  pressure  height. 
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A  second  limitation  to  these  wind  components  is  the  restriction  of 
extremely  large  values  of  dP/dt  that  are  terrain  induced.  After  repeated 
experiments  it  was  determined  that  vertlc  ’  displacements  of  50  mb/3  hr  or 
larger  produce  undesirable  effects.  These  large  vertical  displacements  lead 
to  copious  amounts  of  precipitation  in  areas  of  steep  ascent,  with 
compensating  areas  of  extreme  dryness  on  the  descent  side.  Additionally,  large 
values  of  upward  motion  aggravate  an  inherent  tendency  toward  dryness  in  the 
cloud  model.  This  aridity  is  caused  primarily  by  the  lack  of  low-level 
moisture  sources,  such  as  the  Gulf  of  Mexico,  being  included  in  the  model. 

To  correct  these  def lcleneles ,  the  vertical  displacement  at  each  grid  point 
is  chocked  for  values  greater  than  50  mb/3  hr.  When  this  criteria  is 
exceeded,  each  of  the  horizontal  trajectory  components  is  halved.  A  new 
vertical  displacement  is  computed  from  these  shortened  trajectories,  and  tills 
process  is  continued  until  the  extreme  values  of  DP/dt  are  reduced  below  the 
50  mb/3  hr  criteria. 

3. A. 2  Lateral  Kffects 

The  horizontal  grid  on  which  these  calculations  are  made  Is  superimposed 
on  a  Northern  and  Southern  Hemispheric  polar- s tereographic  map  projection. 

The  orientation  of  the  grid  system  on  the  map  projection  is  such  that  10°K  and 
170°W  longitudes  arc  parallel  to  the  x  axis.  Figures  1  and  2  show  the  two 
hemispheric  projections  and  the  lateral  boundary  of  the  grid  system. 

Generally,  the  area  that  is  inscribed  by  this  octagon  shaped  lateral  boundary 
i;  poleward  of  10°  latitude  excluding  most  of  the  Tropics.  The  model  computes 
trajectories  for  all  grid  points  that  are  Internal  to  this  boundary. 

The  prime  lateral  constraint  on  the  trajectory  computations  is  to  prevent 
origins  from  being  outside  the  octagon  (lateral  boundary).  Generally, 
ail  trajectories  within  two  whole  mesh  (200  tun  between  grid  points)  grid 
points  are  checked  to  Insure  that  this  condition  is  preserved.  Specifically, 
the  actual  constraint  depends  upon  the  trajectory's  proximity  to  a  particular 
boundary.  For  the  grid  rows  and  columns  that  are  adjacent  to  the  boundaries 
that  are  parallel  to  the  x-  and  y  axes,  the  model  does  not  permit  x's  and 
y's  that  are  outside  the  boundary.  Those  parcels  originating  outside  the 
lateral  boundaries  are  truncated  to  the  lateral  boundary. 

Implementation  of  this  truncation  concept  near  diagonal  boundaries  is  more 
complicated.  The  following  mathematical  section  formally  describes  the  system 
used.  Conceptually,  the  diagonal  is  rotated  45°,  the  truncation  concept  Is 
applied  and  then  the  diagonal  13  rotated  back  to  the  initial  orientation. 

Wiiile  the  math  appears  difficult,  the  results  are  straightforward. 

In  order  to  restrict  the  trajectory  origins  that  are  computed  In  the 
it  lnlty  of  the  diagonal  boundaries,  a  modified  x,y-coordinate  system  needs 
L  be  established.  First,  a  unit  vector  perpendicular  to  each  of  the  diagonal 
b  undaries  is  defined.  Using  Figure  9  and  geometric  considerations  the  unit 
vector  u  can  be  specified  in  terms  of  1  and  j  as: 
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(10) 


u  =  '  0.707  L  i  1  0.7071  j 


whore  1  and  J  are  unit  vectors  parallel  to  the  x-  and  y-axes  respectively  and 
the  */-  signs  determine  which  of  the  four  diagonal  boundaries  is  being 
considered .  To  determine  if  the  horizontal  vector  D  has  a  component  that  is 
along  u  and  directed  inward,  the  angle  between  "if  and  "if  is  calculated.  Since, 
by  definition 


u  •  D  =  |uj  | D |  cos  0 


D  =  Axi  +  Ayj 


then,  substituting  for  u  and  D,  the  formula  for  cos  ©  la: 

~  o . 

(x  +  y  ) 

Therefore,  if  u*D  is  greater  than  zero,  then  cos  ©  is  greater  than  zero  and  D 
has  a^component  along  u  and  directed  inward.  If  this  criteria  is  not  met, 
then  D  is  modified  by  subtracting  the  jj^lue  of  the  normal  t^o  the  edge  of  the 
octagon  component  that  is  parallel  to  N.  This  component,  N,  is  defined  as: 

| N {  =  j  D |  cos  0  ^  ^ 

substituting  for  cos  9  from  (11),  N  is: 


Using  (12^  for  D  and  (15)  for  7f,  the  component  of  D  parallel  to  the  diagonal 
boundary  M  is  calculated  from  the  vector  subtraction 

M  =  0  -  N  (16) 

The  horizontal  components  of  the  three-dimentlonal  trajectory  are  set  equal  to 
the  value  of  M.  The  effect  of  these  lateral  constraints  is  to  permit  flow 
parallel  to  the  boundary  and  parcel  origins  that  are  internal  to  the  grid 
system,  while  eliminating  flow  from  outside  the  boundary.  Classically  this 
situation  has  been  described  as  a  free-sllp,  non-permeable  boundary  condition. 
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4.  FORECASTING 


The  forecasting  of  temperature  and  moisture  Is  accomplished  after  the 
Initial  fields  of  these  parameters  and  the  forecast  trajectories  are 
determined.  Despite  the  detail  In  the  following  subsections,  the  general 
methodology  of  a  Lagranglan  advectlon  scheme  Is  simple.  This  scheme  uses 
trajectories  obtained  from  wind  forecasts  to  determine  the  origin  of  an  air 
parcel.  The  Initial  value  of  some  weather  element  associated  with  this  air 
parcel  Is  determined.  Then,  the  Initial  value  Is  modified  as  the  parcel 
traverses  the  path  of  the  trajectory.  The  following  subsections  discuss 
certain  atmospheric  processes  which  Influence  the  temperature  and  moisture 
content  of  an  air  parcel  during  Its  Journey. 

4.1  Moisture 

The  forecasting  techniques  applied  to  CPS  consider  those  atmospheric 
processes  which  greatly  affect  cloud  formation  and  dissipation.  Condensation 
pressure  spread  can  be  modified  by  horizontal  advectlon  as  well  as  vertical 
advectlon.  In  addition,  entrainment  Is  especially  Important  near  the 
boundaries  of  air  masses.  Formulations  have  been  developed  at  AFGWC  to 
diagnose  favorable  entrainment  situations  and  then  to  parameterize  this 
effect.  Another  cloud  producing  mechanism  which  must  be  considered  Is  solar 
heating  of  the  Earth's  surface. 

One  of  the  observable  effects  of  the  heating  cycle  Is  the  diurnal  Increase 
in  cloudiness  over  land  by  day  and  subsequent  decrease  at  night.  Details  of 
this  process  and  Its  simulation  will  be  discussed  later.  Finally,  clouds  are 
Influenced  by  moisture  source  regions  as  they  traverse  open  water.  The 
effects  of  precipitation  and  evaporation  on  cloud  formation  will  be  discussed 
In  a  later  subsection. 

4.1.1  Horizontal  Displacement 

The  numerical  characteristics  of  the  forecasting  routines  can  have  a  largo 
effect  on  the  cloud  forecasts.  Trajectory  origins  are  computed  for  each  grid 
point  at  every  time  step.  The  Interpolation  scheme  used  to  define  the  value 
of  the  weather  elements  to  be  forecast  at  these  origins  Is  very  Important 
because  it  Is  applied  each  time  step.  Any  Inherent  error  in  the  scheme 
contributes  addltlvely  to  forecast  Inaccuracy. 

Selecting  the  Interpolation  scheme  that  introduces  the  least  error  Is  not 
a  simple  task.  Most  schemes  assume  that  a  parameter  varies  continuously 
between  known  values.  Linear  Interpolators  assume  linear  variation.  Bessel 
interpolators  compute  the  gradient  (first  derivative)  and  the  rate  of  gradient 
change  (second  derivative)  to  determine  the  unknown  value.  The  discontinuous 
distribution  of  clouds  In  the  atmosphere  makes  the  application  of 
interpolators  to  cloud  fields  very  difficult. 


Despite  this  difficulty,  a  two  point  and  an  eight  point  linear 
Interpolator  are  used  in  5LAYKR.  Each  of  these  Interpolators  are  shown  in 
Figure  10.  The  two  point  Interpolator  rounds  the  origin  location  to  the 
nearest  grid  point.  Linear  interpolation  is  done  between  the  closest 
pressure  following  levels  as  shown  in  Figure  10A.  The  two  point  Interpolator 
can  be  rapidly  computed;  however,  its  accuracy  is  reduced  because  of  rounding 
the  trajectory  length.  The  eight  point  interpolator  reduces  this  inaccuracy 
by  maintaining  the  proper  trajectory  length  (see  Figure  10B) .  However,  It  is 
a  costly  computer  routine  because  interpolation  is  applied  at  each  grid  point 
(9393  grid  points  per  level),  at  each  level,  and  at  every  time  step. 

The  5LAYKR  model  combines  these  two  Interpolators  in  order  to  derive 
benefits  from  each.  The  two-point  interpolator  is  used  at  a  point  when  the 
cloud  analysis  has  50  percent  or  less  cloud  coverage.  The  eight-point  scheme 
is  used  if  the  analysis  has  greater  than  50  percent  coverage.  The  benefits 
are  a  precise  description  of  large  cloud  coverage  (greater  than  50  percent) 
and  reduce  computational  time.  Combining  these  two  interpolators  causes 
different  cloud  advection  rates.  Scattered  cloud  coverage  (50  percent  or 
less)  will  not  be  adveeted  at  the  same  rate  as  broken  or  overcast  coverage 
(greater  than  50  percent).  This  effect  is  not  considered  a  serious  problem. 
Rounding  trajectory  lengths  over  many  time  steps  reduces  this  error  source 
since  some  trajectories  will  be  lengthened  and  others  shortened. 

4.1.2  Vertical  Displacement. 

The  effect  of  vertical  motion  on  cloud  formation  is  not  thoroughly 
understood.  Unsaturated  rising  air  cools  due  to  adiabatic  expansion  until 
water  vapor  condenses  to  form  clouds.  Conversely,  adiabatic  compression  dries 
air  parcels  until  droplets  disappear.  To  simulate  this  process,  the  vertical 
trajectory  component , &  p ,  is  applied  directly  to  the  CPS  value.  This  is 
described  by 

Cp  -  Ci  4  ^p  (17) 

where  Cp  is  the  forecast  CPS  value  and  Cj  is  the  initial  CPS  value.  The 
value  of  Cj  is  determined  at  the  trajectory  origin  by  three  dimensional 
Interpolation  of  the  Initial  CPS  field.  In  (17)  rising  motion  (negative  p) 
when  applied  to  the  CPS  analysis,  produces  larger  cloud  amounts  and  smaller  CPS 
values.  The  convenience  of  representing  clouds  in  terms  of  units  of  pressure 
is  again  obvious. 

There  are  some  limitations  in  the  application  of  (17).  After  advection 
along  the  trajectory,  Cp  is  limited  to  be  greater  than  or  equal  to  7.ero. 
Negative  values  would  indicate  supersaturation,  in  such  cases  the  excess 
moisture  is  condensed  and  added  to  the  quantitative  precipitation. 


4.1.3  Entrainment 


One  problem  of  true  Lagranglan  advectlon  la  that  there  can  be  no 
Interaction  between  parcels.  In  an  attempt  to  parameterize  this  Interaction, 
a  scheme  was  developed  to  mix  environmental  air  with  the  cloud  mass. 
Entrainment  is  especially  needed  near  the  edges  of  clouds  and  cloud- free 
regions  where  modification  due  to  mixing  is  strong. 

The  model  simulates  two  types  of  mixing.  The  first  is  general  entrainment 
applied  to  all  grid  points  while  the  second  is  a  more  specific  formulation 
that  Is  applied  primarily  at  cloud  edges.  The  general  entrainment  provides 
for  large-scale  mixing  of  moving  cloud  systems.  This  simulation  modifies  the 
CPS  value  that  arrives  at  the  terminal  point  of  the  trajectory  by  applying  the 
CPS  that  existed  at  the  grid  point  In  the  previous  time  step.  The  general 
entrainment  is  given  by 


Cm  --  <3Cp  4  Cp)/4  (18) 

where  Cp  Is  the  CPS  at  the  grid  point  from  the  previous  time  step,  Cp  is 
the  forecast  CPS,  and  CD  is  the  final,  mixed  CPS.  This  technique  typically 
works  well  for  the  advectlon  of  large-scale  cloud  systems  into  arid,  desert 
regions  such  as  the  Sahara  Desert  where,  without  this  technique,  clouds  are 
forecast  more  frequently  than  observed. 

A  disadvantage  of  this  technique  Is  that,  In  time,  the  cloud  boundaries 
tend  to  lose  their  dLstlnct  edges.  To  counteract  this  effect,  a  more  specific 
entrainment  method  was  developed  to  preserve  well-defined  cloud  boundaries. 

The  intent  of  this  procedure  Is  to  emphasize  the  Importance  of  vertical  motion 
to  cloud  formation  or  dissipation;  the  method  has  been  researched  and 
emphasized  by  Edson  (1965).  This  technique  Is  more  limited  because  It  Is 
applied  only  when  descending  parcels  cause  cloud-free  grid  points  to  become 
cloudy  or  when  ascending  parcels  cause  cloudy  points  to  become  clear.  This 
situation  Is  represented  in  Figure  11.  The  methodology  of  this  procedure  Is 
to  first  diagnose  whether  the  specific  conditions  are  occurring  and  second  to 
apply  an  appropriate  "wet"  or  "dry"  weighting  factor.  If  a  cloud- free  grid 
point  becomes  cloudy  while  undergoing  descending  motion,  then  an  area-weighted 
state  of  dryness  is  calculated.  Based  on  the  actual  number  of  horizontally 
adjacent  cloud- free  grid  points  and  the  strength  of  the  descent,  the  drying  of 
the  cloud  parcel  is  forced  to  accelerate.  In  a  similar  manner,  should  a  cloudy 
grid  point  become  clear  while  undergoing  rising  motion,  the  surrounding 
moisture  and  the  upward  motion  are  used  to  accelerate  cloud  formation.  As 
Intended,  this  procedure  emphasizes  the  effect  of  vertical  motion  on  cloud 
formation. 
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Pig.  11.  Situations  for  entrainment  near  cloud  boundaries.  In  A,  the 
grid  point  becomes  cloudy  at  the  time  t  +  ^  t  under  descending  motion. 
In  B,  clearing  Is  associated  with  rising  motion. 


-29- 


I.a^V.^XVZ.T'. 


15'  \*t  1*1 1*!  w  I  . 1 1 1 1  f  1(1 1 IJ  HI  IM  W  M II  lyp.f  H  M  M  W  ■  y  M  »iV  »,» ■■  v  w 

hV 

fc£j 

t£ 


'MVWW 


^  PVJ  IP? 


r-Tv 


K--: 


t.r  • * 

l  ; 


4.1. 4  Diurnal  Effects. 

Another  important  cloud  producing  mechanism  is  the  diurnal  solar  heating 
cycle.  With  the  correct  proportions  of  moisture  and  solar  heating,  air 
parcels  rise  from  the  surface,  cool  adiabatlcally ,  and,  if  the  lifting 
conuensation  level  is  reached,  form  clouds.  This  mechanism  is  especially 
important  during  the  summer  months  when  heating  is  strong  and  horizontal 
adveetlon  is  weak.  Because  of  the  parameterization  procedure  described  in  the 
following  paragraphs,  not  only  is  the  solar  heating  mechanism  Incorporated, 
but  also  any  persistant  localized  effects  such  as  those  caused  by  steep 
terrain  features. 

The  diurnal  cloud  forecasting  technique  involves  maintaining  a  two- day 
average  value  of  CPS  at  each  base  time  (0000  GMT  plus  every  three  hours)  and 
for  every  grid  point.  Therefore,  the  difference  between  two  of  these  fields 
at  a  grid  point  represents  a  diurnal  change  of  CPS.  Since  each  mean  field  nas 
a  specific  base  time,  then  the  diurnal  change  is  valid  over  that  period.  For 
example,  if  the  0000  GMT  field  is  subtracted  from  the  0300  GMT  field,  then  the 
difference  represents  the  diurnal  fluctuation  of  clouds  over  that  time 
period.  Since  these  fluctuations  are  calculated  for  each  grid  point,  local 
effects  are  also  incorporated. 

However,  the  strength  of  the  diurnal  effect  may  be  weak  compared  to  other 
meteorological  effects.  For  example,  adveetlon  of  clouds  may  be  the  major 
factor  that  determines  a  cloud  forecast,  with  diurnal  effects  having  a  minor 
role.  In  order  to  diagnose  the  relative  strength  of  the  diurnal  effect,  each 
point  of  the  current  analysis  is  compared  with  the  equivalent  points  of  the 
corresponding  mean  field.  Based  upon  tl  amount  of  disagreement  between  the 
current  analysis  and  associated  running  mean  field,  a  weighting  factor 
ranging  from  one  to  zi  ro  is  assigned.  If  the  analysis  and  the  mean  field 
are  very  similar  at  a  point,  a  weighting  factor  of  1.0  is  assigned  and  the 
total  diurnal  effect  is  Incorporated.  Smaller  and  smaller  weight  factors  are 
assigned  to  each  point  as  differences  between  the  current  analysis  and  the 
mean  field  increase.  A  zero  weighting  factor  is  assigned  if  the  difference 
between  the  analysis  and  the  mean  field  exceeds  a  previously  established  upper 
limit.  For  CPS,  the  value  of  this  criteria  varies  due  to  the  non-linear 
relationship  between  cloud  and  CPS,  as  shown  in  Figure  7.  For  example, 
changes  in  percent  cloudiness  in  the  middle  ranges  of  cloudiness  represent 
smaller  changes  in  CPS  than  cloudiness  changes  at  either  end  of  its  spectrum. 
In  general,  the  diurnal  parameterization  is  not  applied  if  the  difference 
between  the  mean  and  analyzed  fields  is  greater  than  30  percent  change  in 
cloudiness . 

4.1.5  Evaporation  of  Precipitation 

Lastly,  CPS  can  be  modified  by  simulating  low-level  evaporation  of  liquid 
water  that  falls  from  upper  layers.  Liquid  water  amounts  are  in  the  form  of  a 
quantitative  precipitation  forecast  (QPF).  A  discussion  of  this  QPF 
calculation  is  in  section  4.3.4. 
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The  simulation  of  low-level  evaporation  of  falling  precipitation  Is  based 
on  two  observed  conditions  that  affect  evaporation  rates.  First,  very  dry  air 
produces  an  evaporation  rate  greater  than  moist  air.  Second,  relatively  large 
amounts  of  liquid  water  available  for  evaporation  cause  high  evaporation 
rates.  These  conditions  can  be  quantified  so  as  to  produce  new  values  of 
cloudiness.  If  Qj  represents  the  amount  of  QPP  falling  Into  a  layer  (Inches 
per  3  hours)  and  Is  the  initial  CPS,  then  the  new  CPS,  C,2  la 
calculated  using 

C2  *  Ci  -  (Ci  Qj)  (19) 

The  subtraction  of  the  product  In  (19)  Is  necessary  since  low  values  of  CPS 
represent  high  values  of  cloudiness.  In  addition,  should  Q}  equal  one  inch 
per  three  hours,  then  C2  becomes  equal  to  zero;  this  represents  total 
overcast  conditions.  Therefore,  Q^  is  limited  to  values  no  greater  than 
one.  Should  this  condition  occur,  the  value  of  cloudiness  at  the  grid  point 
Is  set  to  be  overcast,  and  the  remainder  of  the  QPP  la  allowed  to  "fall  out." 
this  remaining  QPP,  Q2,  la  determined  by 

<32  *  «1  -  <Ql  <C1  "  C2)/C!)  (20) 

Comparing  (19)  to  the  two  observable  conditions  stated  above  shows  that 
these  conditions  are  now  specified.  Bquation  (19)  allows  large  changes  In  CPS 
(greater  evaporation)  when  either  large  values  of  C^  (dry  air)  or  large 
values  of  Q}  (large  amounts  of  liquid  water)  are  observed.  Bquation  (20) 
shows  that  the  amount  of  QPP  falling  through  the  layer  Is  proportional  to  the 
percentage  change  In  the  CPS  at  the  point.  This  equation  Is  necessary  In 
order  to  determine  the  amount  of  liquid  water  available  for  evaporation  In 
lower  levels  and  to  determine  the  actual  precipitation  amount  reaching  the 
ground.  Both  of  these  equations  are  applied  to  the  two  lowest  grid  levels, 

850  mb  and  gradient. 

4.2  Temperature 

A  Lagranglan  method  for  advectlon  1s  used  to  forecast  temperature,  as  well 
as  to  forecast  CPS  values.  To  produce  a  forecast,  the  temperature  at  the 
initial  point  of  a  trajectory  Is  determined  by  three-dimensional  interpolation 
In  the  initial  temperature  field.  As  the  trajectory  displaces  the  parcel, 
several  atmospheric  processes  are  simulated  to  modify  this  initial  temperature 
value.  The  dry  adiabatic  process  Is  the  most  important,  at  least  over  periods 
of  two  days  or  less.  Other  Important  effects  on  temperature  forecasts,  but 


secondary,  are  the  effects  of  latent  heat  release,  diurnal  heating,  and 
turbulent  mixing.  Finally,  in  order  to  maintain  vertical  consistency,  the 
temperature  lapse  rate  at  each  point  is  checked  to  insure  that  superadlabatlc 
l,ipse  rates  have  not  been  created.  Section  4.2.5  describes  how  superadlabatlc 
conditions  are  handled. 

4.2.1  Dry  Adiabatic  Process 

The  first  law  of  thermodynamics  for  adiabatic  processes  indicates  that  no 
hojt  is  addel  to  or  taken  away  from  an  air  parcel  during  ascent  or  desent  (a 
reversible  process).  This  formulation,  known  as  Poisson’s  equation,  is  written 
as 

Tp  Tj  (Pp/Pj)  °-286  (21) 

where  Tp  is  the  temperature  at  the  termination  of  the  displacement,  Tj  Is 
the  temperature  at  the  origin  of  the  displacement,  and  Pp  and  Pj  are 
pressure  values  at  the  termination  and  origin  of  the  displacement.  The  value 
of  Pp  is  always  the  pressure  level  for  which  a  forecast  is  being  made,  for 
example,  300  mb.  The  exponent,  0.286,  is  the  ratio  R/Cp  where  R  is  the 
specific  gas  constant  for  dry  air,  and  Cp  is  the  specific  heat  capacity  at 
constant  pressure. 

4.2.2  Pseudo  Adiabatic  Process 

Whenever  an  air  parcel  changes  phase  (melts,  freezes,  evaporates, 
condenses,  or  sublimes),  a  quantity  of  heat  must  be  supplied  to  or  taken  away 
from  the  parcel.  This  heat  warms  or  cools  environmental  air  at  a  rate  that 
differs  from  the  dry  adiabatic  (reversible)  process.  To  simulate  this 
pseudo  adiabatic  effect  on  the  temperature  of  saturated  parcels,  a  technique 
using  an  approximate  solution  has  been  developed.  The  ultimate  goal  of  this 
technique  is  to  find  the  temperature  of  a  parcel  that  terminates  at  the  final 
forecast  level,  Pp,  after  it  has  gone  through  a  moist  expansion  process. 

The  approximate  solution  assumes  linear  rates  of  change  of  wet  bulb 
potential  temperature  with  pressure.  This  eliminates  the  need  for 
exponentiation  (an  expensive  computer  operation)  and  produces  accurate  results 
for  the  relatively  small  vertical  displacements  typical  of  three-hour  time 
s  teps . 


4-2.3  Entrainment. 


As  stated  in  Section  4.1.3,  there  needs  to  be  some  mechanism  to  allow 
parcels  of  air  to  interact.  For  temperature,  this  is  especially  critical  for 
two  situations  where  thermal  mixing  is  a  significant  effect  One  involves 
mixing  over  water,  and  the  other  low  level,  turbulent  mixing 

The  importance  of  temperature  modification  of  air  masses  as  they  move  over 
large  water  sources  has  been  well  established.  A  study  by  Reap  (1971)  reports 
the  importance  of  this  modification  and  describes  an  empirical  expression  that 
forecasts  air  temperatures  over  water  regions.  To  produce  an  air  temperature 
forecast,  this  expression  combines  three  fourths  of  the  final  sea- surface 
temperature  plus  one- fourth  of  the  initial  air  temperature.  After 
experimenting  with  several  proportions  at  AFGWC,  the  following 
parameterization  for  air  sea  temperature  interaction  was  selected 

Tm  -  (3Tp  f  TP)/4  (??.) 

where  Tp  is  the  temperature  at  the  grid  point  from  the  previous  time  step, 

Tp  is  the  forecast  temperature,  and  Tpj  is  the  final,  mixed  temperature. 

This  expression  is  applied  at  all  levels  for  all  grid  points  over  water. 

Another  temperature  modification  simulates  the  effects  of  low  level, 
turbulent  mixing.  The  intent  of  this  procedure  is  to  diagnose  conditions 
favorable  to  this  mixing  and  then  apply  the  mixed  condition.  One  such 
condition  that  13  favorable  for  this  mixing  is  parcels  descending  to  the 
gradient  level.  Instead  of  computing  an  adiabatic  change  to  this  level,  the 
standard  atmospheric  lapse  rate  for  the  troposphere  is  used  (5.9  K/100  mb) 

This  lapse  rate  is  less  than  the  adiabatic  one.  Therefore,  parcels  descending 
to  the  gradient  level  are  cooler  than  if  they  descend  at  the  adiabatic  rate. 
The  selection  of  the  standard  atmospheric  lapse  rate  13  solely  based  on  its 
ability  to  produce  gradient- level  temperatures  that  are  cooler  than  those 
produced  by  an  adiabatic  lapse  rate. 

4.2.4  Diurnal. 

Temperature  forecasts  are  also  modified  by  local  diurnal  effects.  The 
method  to  produce  these  effects  is  essentially  the  same  method  as  that  used  to 
dlurnally  modify  CPS  values  (recall  Section  4.1.4).  Two  minor  variations  of 
this  method  need  to  be  noted.  First,  the  mean  temperature  field  is  maintained 
for  eight  days  Instead  of  two.  Second,  the  upper  limit  that  determines 
whether  any  diurnal  effects  are  to  be  applied  is  set  to  283  K.  These  two 
values  have  been  established  after  an  extensive  statistical  evaluation  of 
temperature  forecasts. 


4 .  2 .  5  Superadiabatlc  Lapse  Rate  Modification. 

forecast  temperatures  are  altered  when  the  adveetlon  process  creates 
superadiabatlc  lapse  rates.  These  lapse  rates  are  eliminated  by  adjusting  the 
forecast  lapse  of  temperature.  Klimlnatlng  superadlabal Ic  lapse  rates  Is 
justified  since  they  rarely  would  be  observed  and  only  for  short  durations. 

The  general  procedure  to  remove  a  superadiabatlc  layer  is  to  determine  the 
lapse  rate  between  two  pressure  levels.  If  this  rate  is  superadiabatlc,  then 
the  upper  level  Is  warmed  and  the  lower  level  is  cooled.  This  warming  and 
cooling  is  done  in  small  increments,  and  the  sounding  is  continuously  checked 
to  ensure  that  no  superadiabatlc  lapse  rates  are  created  at  other  levels  by 
this  adjustment  procedure. 

4.3  Derived  Quantities 

Temperature  and  CPS  are  the  only  meteorological  elements  that  are  advected 
by  this  forecasting  scheme.  However,  additional  elements  may  be  derived  from 
forecasts  of  them.  These  non  advected,  secondary  forecast  fields  Include 
layered  and  total  cloud  cover,  dew  point  depression,  stability,  quantitative 
precipitation  forecasts,  cloud  type,  and  icing.  The  following  sections 
describe  the  methodology  that  is  used  to  obtain  these  derived  fields 

4.3.1  Layered  and  Total  Cloud  Cover 

Layered  cloud  forecasts  are  derived  from  forecast  CPS  fields  using  tables 
(Appendix  B)  that  convert  CPS  to  cloud  amounts.  P’rom  these  layered  cloud 
forecasts,  a  total  cloud  field  is  created  by  the  technique  discussed  In 
Section  2.2.2. 

The  choice  of  CPS  as  the  adveetlon  element  for  cloud  forecasting  Is  not  a 
simple  one.  For  example,  primitive  equation  forecast  models  that  incorporate 
moisture  can  Initialize  and  forecast  a  moisture  related  element  such  as 
specific  humidity  or  dew  point  depression  (T  T^).  Areas  of  clouds  are  then 
Inferred  from  these  moisture  fields.  In  fact,  early  cloud  forecasting  models 
at  AFGWC  used  dew  point  depression  to  infer  cloud  amounts.  To  determine  which 
meteorological  element  produces  the  most  realistic  cloud  forecast,  an 
experiment  was  conducted  at  AFGWC  during  the  summer  of  1974  contrasting  CPS 
and  dew  point  depression. 

fiie  experiment  consisted  of  independently  initializing,  forecasting,  and 
verifying  dew  point  depressions  and  CPS.  An  automated  verification  program 
compared  the  following  criteria:  forecasts  and  observations  of  dew  point 
depress  ions ,  forecasts  of  CPS  (converted  to  dew  point  depression)  to 
observations  of  T  Td,  forecasts  of  CPS  to  observed  clouds,  and  forecasts  of 
dew  point  depression  to  observed  clouds.  While  the  statistical  results  of 


this  experiment  are  not  available,  the  qualitative  results  were  summarized  as 
follows:  In  verifying  dew-point  depression  forecasts  against  station 
observations  of  dew-point  depressions,  advected  MULTAN  dew-point  depressions 
(a  scheme  that  initializes  and  forecasts  dew-point  depressions)  shows  more 
skill  almost  everywhere  than  CPS  (a  scheme  that  Initializes  and  forecasts 
CPS).  However,  advected  MULTAN  dew-point  depression  forecasts  are  poorly 
correlated  with  cloudiness.  Because  of  these  results  and  the  operational 
requirement  for  a  cloud  forecasting  model,  CPS  was  chosen  to  be  the  advectable 
moisture  element  in  this  forecasting  scheme.  Since  CPS  is  the  only  advectable 
quantity,  the  verification  statistics  should  correlate  well  with  cloud 
observations  but  not  with  dew-point  depressions.  Bven  though  forecasts  of 
dew-point  depression  are  produced,  they  are,  in  fact,  derived  from  forecasts 
of  CPS  using  the  method  described  in  the  following  section.  Users  of  APOWC 
dew-point  depression  forecasts  should  be  aware  of  the  procedures  that  derive 
these  forecasts. 

4.3.2  Dew-Point  Depression 

Forecasts  of  dew-point  depression  (DPD),  T-Td,  are  produced  directly 
from  the  forecast  CPS  values.  The  equation  for  the  conversion  from  CPS  to  DPD 
is : 

DPD  --  CPS  /  [A0(P)  +  Al  (P)  «  CPS  )  (23) 

Aq(P)  =■  Aqq  ♦  Aqi  *  P 
Ai(P)  =  A10  ♦  An  •  P 

where 


A00  =  1.46917 
Aq!  =  1.363058-2 

and 


A10  -  -9.011778-3 
An  =  1.7772  8-6 

The  resulting  error  is  less  than  5%,  when  compared  with  an  exact  conversion. 
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4.3.3  Stability 


Using  temperature  and  moisture  forecasts,  a  measure  of  the  atmosphere's 
static  stability  can  be  determined.  For  those  grid  points  that  have  a  terrain 
height  below  the  850  mb  level,  a  standard  Showalter  Stability  Index  (Huscke, 
1959)  is  computed.  Using  the  850  mb  temperature  and  dew  point,  the  parcel's 
pressure  and  temperature  at  its  lifted  condensation  level  (LCL)  is 
determined.  The  parcel  at  the  LCL  is  lifted  moist  adiabatically  to  500  mb 
where  the  parcel’s  temperature  is  compared  to  the  500  mb  environmental 
temperature.  The  difference  between  these  two  temperatures  is  the  Showalter 
Stability  Index.  The  moist  adiabatic  computation  is  the  one  mentioned  in 
Section  4.2.2. 

When  the  terrain  is  above  850  mb,  a  different,  more  empirical  approach  is 
taken.  Based  upon  observational  evidence,  a  technique  that  produces  a 
continuous  field  of  stability  indexes  in  a  steep  terrain  region  is  used.  It 
was  noted  that  the  difference  between  the  700  mb  and  500  mb  wet  bulb  potential 
temperature  is  a  stability  index  that  is  essentially  the  same  as  if  a  standard 
Showalter  Tndex  was  computed  from  850  mb.  For  those  grid  points  which  have 
terrain  heights  above  700  mb,  no  stability  index  is  forecast. 

4.3.4  Quantitative  Precipitation  Forecasts 

The  technique  to  calculate  a  quantitative  precipitation  forecast  (QPF)  is 
similar  to  the  formulation  that  is  reported  by  Haltiner  (1971).  This  modified 
formula  is 


Q 


t+1 


dz 


dc 


(24) 


where  Q  is  QPF,  w  is  the  saturation  mixing  ratio  of  the  parcel,  and  e  is  the 
standard  density  for  the  layer.  The  Integration  is  over  time  step, At,  and 
over  the  thickness  of  a  layer  zu  z^.  The  vertical  summation  of  (24)  is 
done  independently  at  the  500  mb,  700  mb,  and  850  mb  levels  The  upper  level, 
zu,  refers  to  a  point  midway  between  the  standard  pressure  level  being 
forecast  and  the  next  standard  level  above  it.  The  same  condition  is  true  for 
zj,  except  it  is  below  the  forecast  level.  The  amount  of  precipitation  that 
reaches  the  surface  is  the  sum  of  the  forecasts  for  the  independent  layers. 


Generally,  the  calculation  of  each  of  the  terms  in  (24)  consists  of 
manipulating  quantities  that  either  have  been  determined  or  using  the  methods 
discussed  above.  The  first  variable  is  the  change  in  the  saturation  mixing 
ratio,  dw,  which  is  the  difference  between  the  mixing  ratio  at  the  start  and 
the  end  of  the  moist  adiabatic  expansion.  The  difference  represents  the 
amount  of  liquid  water  that  condenses  and  is  expressed  in  units  of 
gm  H20/kg  air.  All  liquid  water  is  assumed  to  fall  as  precipitation.  Using 
equations  derived  from  the  definition  of  the  saturation  mixing  ratio  (w 
we),  the  mixing  ratio  at  the  start  (subscript  s)  and  at  the  end  (subscript 
e)  of  the  moist  expansion  is 


wg  =  0.622 

( e3/Pg ) 

(25) 

we  =  0.622 

( ee/P e) . 

(26) 

The  term  P3  is  the  pressure  of  the  LCL  and  Pe  is  the  standard  pressure 
level  for  which  the  forecast  is  being  made. 

Saturation  vapor  pressure  (es  and  ee  in  the  above  equations)  is 
calculated  in  the  manner  described  by  Lowe  and  Ficke  (1974).  Their  technique 
involves  solving  a  sixth  -  order  polynomial  approximation  that  specifies  e  as 
a  function  of  temperature  only.  The  temperatures  for  this  polynomial  are 
those  at  the  beginning  and  end  of  the  moist  expansion.  The  advantage  of  using 
a  polynomial  solution,  rather  than  one  requiring  exponentiation,  is  a  faster 
computation  of  the  variable.  As  discussed  by  Lowe  and  Ficke,  accuracy  of  the 
approximation  is  very  good. 

To  this  point  of  the  discussion,  it  has  been  assumed  that  total  saturation 
occurs  exactly  at  the  LCL.  However,  condensation  and  subsequent  precipitation 
have  been  shown  to  occur  at  relative  humidity  values  les3  than  100  percent 
when  averaged  over  a  large  area.  To  allow  for  this  observation,  the  terms  in 
(24)  are  not  evaluated  at  the  LCL,  but  rather  at  a  level  below  it.  It  has 
been  determined  experimentally  at  AFGWC  that  this  level  is  7  mb  below  the 
LCL.  This  level  i3  identical  to  a  CPS  value  of  7  mb,  and  it  corresponds  to  a 
cloudiness  of  approximately  95  percent. 

The  remaining  terms  in  (24),  'e'  dz  ,  together  represent  the  weight  of  air 
per  unit  area.  When  these  terms  are  combined  with  w,  the  desired  units  of 
weight  of  water  per  unit  area  are  obtained.  The  term,  z,  is  computed  using  a 
form  of  the  hydrostatic  equation. 


z  =  R  Te  In  (Pe/Pu)/g 


(27) 
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where  R  is  the  gas  constant  and  g  is  the  acceleration  due  to  gravity.  As 
reported  by  Haltiner  (1971),  the  precipitation  equation  (24)  gives  good 
results  for  synoptic- scale  pressure  systems  in  mid-latitudes  but  it  is 
inadequate  in  areas  of  convection,  especially  in  tropical  systems. 

4-3-5  Cloud  Type 

Techniques  have  been  developed  to  combine  temperature  and  CPS  in  order  to 
classify  the  types  of  cloud  structures.  The  primary  distinction  between  cloud 
types  is  the  distinction  between  stratiform  and  cumuliform  clouds.  A 
stability  value  is  calculated  for  the  850  mb,  700  mb,  and  500  mb  pressure 
levels  to  assist  in  this  distinction.  This  stability  determination  consists 
of  comparing  wet  bulb  potential  temperatures,  between  two  levels.  This 
procedure  is  identical  to  the  description  in  Section  4-3.2.  For  example,  if 
the  0W  at  850  mb  is  greater  than  the  9W  at  700  mb,  a  positive,  upward 
buoyancy  exists  at  the  700  mb  level.  Should  clouds  be  forecast  at  a  pressure 
level,  positive  buoyancy  at  that  level  implies  cumuliform  clouds,  whereas 
negative  buoyancy  implies  stratiform.  Vertical  cloud  structures  such  as 
towering  cumulus  and  cumulonimbus  are  determined  by  applying  these  checks  at 
all  levels.  For  example,  if  the  850  mb  level  is  unstable  and  forecast  clouds 
exist,  the  700  mb  cloud  and  stability  are  checked.  If  clouds  and  upward 
buoyancy  exist  at  700  mb,  then  the  500  mb  level  is  checked.  Should  upward 
buoyancy  exist  at  500  mb,  the  cloud  is  typed  as  a  cumulonimbus,  but  if 
downward  buoyancy  exists,  the  cloud  is  classified  as  towering  cumulus.  Some 
typical  combinations  of  stability  and  clouds  are  shown  Figure  12. 

4.3.6  Precipitation  Type 

The  primary  reason  for  classifying  the  types  of  cloud  structures  is  to 
provide  a  basis  for  determining  precipitation  type.  Forecasts  of 
precipitation  types  are  categorized  as  either  continuous  (from  stratiform 
clouds)  or  showery  (associated  with  cumuliform  clouds).  The  development  of 
criteria  to  determine  precipitation  types  and  their  intensities  was  done  by 
comparing  forecast  values  of  cloud  type,  cloud  cover,  and  static  stability  to 
precipitation  observations.  Results  of  these  comparisons  indicate  that  no 
precipitation  is  observed  when  the  total  cloud  cover  is  less  than  3/8  (37 
percent)  coverage  at  a  grid  point.  It  should  be  remembered  that  values  at 
grid  points  represent  average  conditions  over  a  large  horizontal  area. 
Therefore,  these  precipitation  forecasts  should  not  be  viewed  as  a  point 
forecast . 

For  showery  precipitation  (rainshower,  snowshower,  and  thunderstorm)  a 
cloud  type  of  towering  cumulus  or  cumulonimbus  is  required.  Also,  a  minimum 
cloud  cover  of  25  percent  is  necessary  at  either  the  850  or  700  mb  levels. 

The  distinction  between  rainshowers  and  thunderstorms  is  based  primarily  upon 
whether  the  cloud  type  is  towering  cumulus  or  cumulonimbus.  Snow  showers 
occur  if  the  gradient  level  temperature  is  less  than  273°  K. 
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Fig.  12.  Some  typical  cloud  types  associated  with  various  buoyancies. 
Buoyancies  are  determined  by  the  layer  difference  in  wet-bulb  potential 
temperature.  Unstable  (  +  )  conditions  exist  if  the  wet- bulb  potential 
temperature  at  the  lower  layer  is  greater  than  that  at  the  higher 
layer.  Stable  (  )  conditions  exist  if  the  reverse  is  true.  Due  to 
having  information  only  at  discrete  layers,  it  is  difficult  to  make 
cloud  top  or  base  estimates. 


The  Intensity  of  rainshowers  is  inferred  from  the  700  mb  static  stability, 
whereas  ralnshower  intensity  associated  with  thunderstorms,  is  determined  by 
computing  a  total-totals  stability  index.  This  index  is  defined  by  Miller 
(1972).  The  threshold  value  differentiating  moderate  or  heavy  rainshowers 
associated  with  thunderstorms  has  been  correlated  with  observational 
evidence.  The  reader  should  note  that  shower-producing  mechanisms  such  as 
steep  terrain  and  diurnal  heating  are  not  considered. 

For  continuous  precipitation  (rain,  drizzle,  and  snow)  a  large  areal 
coverage  of  cloudiness  is  required.  Cloud  types  other  than  towering  cumulus 
or  cumulonimbus  and  cloud  cover  greater  than  75  percent  at  850  mb  or  700  mb 
are  the  minimum  criteria  for  this  precipitation  type.  Snow  is  differentiated 
from  liquid  precipitation  whenever  the  gradient  level  temperature  is  less  than 
273°  K. 

The  intensity  for  continuous  precipitation  is  a  function  of  the  moisture 
depth  and  the  buoyancy  at  700  mb.  Por  example,  heavy  continuous  rain  is 
forecast  when  cloud  amounts  are  greater  than  75  percent  at  all  three  levels 
(850,  700,  and  500  mb)  and  the  850  mb  wet-bulb  potential  temperature  is  more 
than  3  K  greater  than  the  700  mb  value.  Lesser  intensities  are  forecast  for 
shallower  cloud  depths  and  smaller  buoyancy  values. 

4.3.7  Icing 

Forecasts  of  icing  can  be  derived  using  temperature  and  moisture 
information.  The  criteria  to  forecast  icing  type  and  intensity  are  the  same 
as  that  in  AVTS/TR  80/001,  Forecaster's  Guide  on  Aircraft  Icing.  Table  1  is  a 
summary  of  these  criteria.  Initially,  all  icing  intensities  are  light.  When 
cold  air  advection  greater  than  2°K/3  hr  is  forecast,  the  intensity  is 
increased  to  moderate.  Warm  air  advection  of  0.1°K/3  hr  will  reduce  the 
intensity  to  trace.  For  temperatures  outside  of  the  ranges  in  Table  1,  no 
icing  is  forecast.  All  icing  intensities  are  Initially  classified  as  light. 


Table  1.  Icing  Forecast  Criteria. 


TEMPERATURE  (K) 

DEW  POINT  DEPRESSION  (K) 

ICING  TYPB 

CLOUD  TYPB 

272-266 

-2 

Clear 

Cumulus 

Rime 

Stratus 

265  258 

-3 

Mixea 

Cumulus 

Rime 

Stratus 

257  251 

-4 

Rime 

All  Clouds 
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5.  HIGH  RESOLUTION  CLOUD  PROGNOSIS 


5.1  Model  Specifications 

The  High-Resolution  Cloud  Prognosis  (HRCP)  module  forecasts  percent 
cloudiness  for  15  levels  on  an  eighth-mesh  grid.  This  eighth-mesh  grid 
coincides  with  the  one  used  in  the  RTNEPH  (Fig.  13).  In  addition,  a  percent 
total  cloudiness  is  forecast  for  each  eighth-mesh  point.  Forecasts  are  made 
in  three-hour  time  steps  to  a  maximum  of  9  hours. 

The  HRCP  is  capable  of  making  a  forecast  for  that  portion  of  the 
eighth-mesh  grid  within  the  AFGWC  octagon.  However,  the  advection  module  of 
the  HRCP  only  forecasts  for  a  single  window  of  this  area  in  any  one  computer 
run.  This  window  can  be  located  anywhere  in  the  octagon;  however,  the  window 
must  coincide  with  one  or  more  of  the  boxes  of  the  eighth- mesh  grid.  These 
boxes  are  numbered  from  one  to  64  as  shown  in  Fig  13. 

5.2  Options 

The  HRCP  can  be  initiated  for  any  starting  hour,  and  it  phases  into  the 
5 LAYER  database  to  use  the  trajectories  already  stored  there.  It  has  the 
option  to  make  any  length  forecast  out  to  9  hours  in  three-hour  increments. 

5.3  Forecast  Mechanism 

The  HRCP  is  initialized  directly  from  the  RTNEPH  converted  to  look  like 
the  old  3DNBPH  database  (primarily  15  layers  of  cloud  information).  No 
modification  of  the  cloud  amounts  is  made.  To  produce  an  eighth-mesh 
cloudiness  forecast,  HRCP  uses  information  available  in  the  5LAYBR  and  the 
RTNEPH  cloud  models.  The  first  major  step  is  to  calculate  the  local  change  of 
CPS  with  time  at  the  half-mesh  grid  points.  Trajectories  are  available  only 
on  the  half-mesh  grid,  the  advection  of  the  initial  CPS  field  is  done  on  this 
grid  scale.  The  local  change  of  CPS  results  from  taking  the  difference 
between  the  initial  and  forecast  CPS  fields. 

In  order  to  apply  this  half-mesh,  5LAYBR,  local  tendency  of  CPS,  an 
eighth-mesh,  15-level  format  must  be  used.  To  do  this,  the  CPS  tendency  is 
horizontally  interpolated  to  the  eighth- mesh  grid  points.  This  tendency  is 
not  vertically  interpolated,  but  rather  each  of  the  forecast  five  levels  are 
directly  applied  to  their  corresponding  levels  of  the  15-level  grid  (cloud 
analysis).  For  example,  the  tendency  at  300  mb  is  applied  to  layers  13,  14, 
and  15  (Fig.  14). 

A  slightly  different  methodology  is  used  for  the  gradient  level.  The  CPS 
tendency  for  this  level  is  applied  only  to  levels  five  and  six.  The  cloud 
amounts  in  levels  one  to  four  are  persisted.  The  intent  of  this  method  is  to 
eliminate  the  advection  of  very  low-level  clouds.  Advection  of  low  clouds 
such  as  fog  would  be  a  questionable  procedure. 

Lastly,  this  CPS  tendency  is  applied  directly  to  the  3DNEPH  to  produce  an 
eighth-mesh,  15-level  cloud  forecast.  This  process  is  depicted  in  Fig.  15. 

In  addition,  a  total  cloud  forecast  is  produced  by  stacking  these  15  levels. 
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LAYER  LEVEL 


Fig.  14.  Association  of  the  15  fixed  levels  (left)  in  HRCP  that  the 
trajectories  in  the  five  forecast  layers  of  5LAYKR  are  applied  to. 
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PROCESSING  PLOW  CHART 


5 LAYER  Forecast  Field 
(1/2  mesh,  5  layers) 


SLAYER  Initial  Eield 
(1/2  mesh,  5  lavers) 


,Cnanae  Field  (CPS) 
1(1/2  mesh.  5  layers) 


Chanqe  Field  (CPS) 
(1/2  mesh,  5  layers) 


Interpolate 


► 


Chanqe  Field  (CPS) 

(1/8  mesh,  15  layers) 


Chanqe  Field  (CPS) 
(1/8  mesh,  15  layers] 

+ 

Initial  3DNEPH  (CPS) 
(1/8  mesh,  15  lavers) 

ss 

HRCP  Forecast  (CPS) 
(1/8  mesh,  15  lavers) 

i 


Convert 


HRCP  Forecast  Cloud 
(1/8  mesh,  15  layers) 


Fig.  15.  Flow  chart  of  HRCP  processing  of  5LAYBR  and  RTNEPH  data 
to  produce  an  eighth  mesh  forecast. 
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6.  TROPICAL  CLOUD  POR8CASTS 


6.1  Model  Specification 

The  Tropical  Cloud  Forecasting  Model  (TRONEW)  produces  forecaste  for  throe 
layers  (low,  middle,  and  high)  on  the  APGWC  Northern  Hemisphere  and  Southern 
Hemisphere  half-mesh  supergrids  from  approximately  25°N  to  25°S.  Forecasts 
are  made  every  three  hours  In  three-hour  increments  out  to  21  hours. 

6.2  Forecast  Mechanism 

TRONBW  uses  the  premise  that  In  the  tropics  there  are  diurnal  fluctuations 
of  cloudiness;  and  therefore,  clouds  that  were  observed  yesterday  at  a  certain 
time  of  day  will  again  be  observed  today  at  that  same  time.  This  premise 
works  well  except  In  the  vicinity  of  moving  tropical  disturbances  (e.g. 
typhoons,  and  easterly  waves). 

To  Initialize  the  half- mesh  grid,  a  9-point  weighted  average  (interior 
grid  points  In  Pig.  4)  of  the  eighth-mesh  RTNBPH  grid  is  used.  This  compacted 
analysis  is  then  saved,  so  that  if  a  six-hour  forecast  is  needed  from  a 
23/0600Z  base  time,  TRONEW  will  use  the  22/1200Z  analysis  as  that  forecast. 


7.  MODEL  EVALUATION 


7.1  Verification 

The  primary  purpose  In  verifying  any  forecast  scheme  Is  to  provide  a  basis 
for  future  model  Improvement.  Improvement  may  occur  by  developing  better 
forecast  techniques  or  by  Identifying  characteristic  model  strengths  and 
weaknesses.  As  new  techniques  are  developed  and  evaluated,  they  are  compared 
to  objective  verification  statistics.  Any  model  degradation  or  Improvement  Is 
Immediately  Identified.  Forecasters  who  use  numerical  guidance  routinely  make 
subjective  evaluations  of  the  model  output.  By  Identifying  situations  that 
typically  produce  good  or  bad  forecasts,  the  forecaster  can  modify  the  output 
and  therefore  improve  the  forecast. 

The  5LAYBR  model  Is  evaluated  by  both  subjective  and  objective 
verification  schemes.  Forecasters  at  AFGWC  routinely  use  5LAYBR  forecasts  to 
assist  them  In  providing  worldwide  meteorological  support.  The  Special 
Support  Division  and  the  Forecasting  Services  Division  are  the  primary  users 
of  the  5LAYBR  model  at  AFGWC.  Special  Support  forecasters  use  5LAYBR  for 
their  Initial  forecast  guidance  and  for  quality  control  of  the  product. 

Quality  control  Is  maintained  by  comparing  verification  of  the  forecaster’s 
product  to  that  of  the  model. 

Forecasting  Services  Division  forecasters  use  5LAYKR  In  a  similar  manner. 
The  Horizontal  Weather  Depiction  (HWD)  Section  of  this  division  Is  responsible 
for  forecasting  clouds.  Icing,  and  precipitation  type  on  a  hemispheric 
coverage.  These  forecasts  provide  the  primary  guidance  for  similar  forecasts 
that  have  more  specific  application.  Por  example,  they  are  used  for  Strategic 
Air  Command  (SAC)  low- level  route  forecasts.  Tactical  Air  Command  (TAC)  range 
forecasts,  and  U.S.,  European,  and  Aslan  terminal  forecasts.  Because  so  many 
forecasts  rely  on  the  "lead"  forecast  of  the  HWD  section,  the  product  Is 
continuously  quality  controlled. 

Objective  verification  of  5LAYKR  forecasts  is  accomplished  by  automated, 
model  verification  programs  that  produce  statistical  evaluation  of  layered  and 
total  cloud  forecasts.  Persistence  "forecasts"  are  also  verified  as  a 
standard  to  compare  model  performance. 

Since  the  RTNRPH  Is  the  standard  against  which  cloud  forecasts  are 
verified,  a  discussion  of  the  model’s  characteristics  is  necessary.  The 
RTNEPH  uses  data  from  all  available  sources.  Satellite,  radiosonde  reports, 
standard  surface  observations,  and  aircraft  reports  are  merged  to  produce  a 
hemispheric  cloud  analysis.  Each  of  these  data  sources  contribute  their  own 
characteristics  in  determining  the  "true”  cloud  amount.  Satellite  data  for 
the  RTNKPH  13  provided  by  the  Defense  Meteorological  Satellite  Program 
( DMSP ) .  This  program  normally  consists  of  two  polar-orbiting  satellites. 
Sensors  on  the  satellite  accurately  determine  the  horizontal  cloud  coverage. 
However,  the  vertical  distribution  of  clouds  must  be  inferred  from  other  data 
sources.  Vertical  temperature  profiles  used  with  infrared  satellite  sensors 
can  give  approximate  vertical  positioning.  Aircraft  reports  and  surface 


observations  also  aide  in  this  determination.  DMSP  satellite  data,  because  of 
polar  type  orbits,  do  not  give  hourly  global  coverage.  Therefore,  at  a 
particular  verification  time  new  satellite  data  may  not  be  available.  Since 
the  RTNBPH  uses  persistence  in  the  absence  of  new  data,  a  forecast  of 
persistence  would  verify  as  being  correct. 

Most  grid  points  in  the  model  verification  program  are  located  over  land 
masses  to  partially  alleviate  any  bias  favorable  to  the  persistence  forecast. 
However,  surface  observations  of  clouds  contribute  their  own  source  of  error. 
The  line-of-sight  of  a  surface-based  observer  is  only  10  to  20  miles  depending 
on  obstructions  and  viewing  angle  with  respect  to  the  terrain.  These 
restrictions  typically  cause  surface  observers  to  overestimate  horizontal 
cloud  distribution.  Low  overcast  cloud  layers  bias  surface  reports  by 
preventing  reports  of  cloud  conditions  above  the  overcast.  The  RTNBPH  persists 
the  cloud  analysis  in  the  absence  of  new  layered  cloud  observations.  As  was 
discussed  for  satellites,  there  continues  to  be  a  favorable  bias  in  the 
verification  toward  the  persistence  forecast.  Despite  these  problems  the 
5LAYBR  model  out  performs  the  persistence  forecast.  A  more  complete 
discussion  of  the  comparative  cloud  forecast  skill  of  5LAYBR  was  prepared  by 
Mitchell  (1982).  Figure  16  provides  a  graphical  representation  of  the  5LAYBR 
performance  during  May  1985  to  April  1986.  The  verification  variable  is  the 
percentage  of  forecasts  of  total  cloud  that  verified  within  +25%  of  the  total 
cloud  cover  analyzed  by  RTNBPH.  As  with  other  forecast  models,  forecast  skill 
deteriorates  with  increasing  length  of  forecast. 

7.2  Five  Layer  Model 

Numerical  models  provide  an  objective  means  of  producing  forecasts.  When 
used  as  a  forecaster’s  aid,  they  serve  as  a  basis  from  which  a  forecaster  can 
produce  an  improved  forecast.  By  learning  inherent  characteristics  of  a 
model,  a  forecaster  can  better  manage  his  use  of  these  aids.  Ideally,  the 
forecaster  should  spend  his  greatest  efforts  on  situations  of  known  model 
deficiencies.  Conversely,  little  time  should  be  spent  on  situations  where  the 
model  performs  well.  The  purpose  of  thi3  subsection  is  to  identify  strengths 
and  weaknesses  of  the  5LAYER  model. 

The  5LAYBR  model  depends  on  other  models  to  provide  cloud  and  temperature 
analyses  and  wind  forecasts.  Characteristics  of  these  models  inherently 
become  a  part  of  the  5LAYBR  model.  Wind  forecasts  from  the  AWS  GSM  model  are 
formed  into  trajectories  as  described  in  Section  3  and  Appendix  D. 
Characteristics  that  are  typical  of  the  AWS  GSM  are  also  typical  of  the 
5LAYBR. 

The  5LAYER  model  produces  the  best  cloud  forecasts  at  mid- tropospheric 
levels  (700  and  500  mb)  and  at  middle  latitudes.  This  should  be  expected 
3 ince  dynamic  wind  models  (like  AWS  GSM)  produce  their  best  forecasts  in  his 
region.  Also,  clouds  that  are  associated  with  moderate  or  strong 
synoptic-scale  weather  patterns  are  identified  and  forecast  with  the  greatest 
skill .  Again,  this  is  to  be  expected  since  mid- latitude  storm  systems  are  of 
a  sufficient  scale  to  be  resolvable  by  the  5 LAYER  grid  system.  These  cloud 
systems  also  have  well  defined  cloud  boundaries  which  are  necessary  for 
accurate  cloud  mapping  and  for  predicting  movement  of  these  boundaries. 
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Fig.  16.  One  year  summary  of  monthly  Northern  Hemisphere  5LAY8R 
total  cloud  25/25  skill  scores  (percentage  of  forecasts  within 
+25%  of  observed  total  cloud)  for  6,  12,  24,  36,  and  48  hour 
forecasts.  Decreasing  skill  with  forecast  length  is  apparent. 
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RTN8PH  data  must  be  put  In  a  workable  format  for  the  5LAYER  model. 

Section  2.2.1  discusses  the  methodology  of  locating  the  vertical  the  cloud 
data  and  representing  25  eighth  mesh  grid  points  by  one  half  mesh  point 
Figure  17  shows  how  horizontal  data  compaction  affects  cloud  prediction  on  the 
larger  scale.  Figure  17  shows  a  cloud  boundary  that  is  averaged  to  a 
half- mesh  grid  point.  Any  averaging  done  across  a  cloud  edge  causes  a 
smoothing  of  that  edge.  Rather  than  a  distinct  change  from  100  percent  to 
zero  percent  cloud  cover,  the  half -mesh  representation  becomes  100  percent,  50 
percent,  and  zero  percent.  For  large  scale  cloud  systems,  this  reduction  in 
cloud  mass  is  not  a  significant  problem.  However,  for  systems  represented  by 
only  a  few  grid  points,  a  large  percent  of  the  cloud  mass  may  be  lost.  This 
los3  is  solely  caused  by  compacting  the  cloud  data. 


Restrictions  imposed  on  the  trajectory  computation  are  discussed  in 
Section  3.  The  effect  of  these  restrictions  on  cloud  prediction  has  been 
observed  and  are  discussed  here.  Trajectories  are  modified  at  the  lower  and 
lateral  boundaries  of  the  5LAYBR  grid  system.  Trajectories  near  lateral 
boundaries  are  restricted  to  be  parallel  to  the  octagon  edge  This 
restriction  causes  clouds  to  be  advected  around  the  grid  edge.  This  procedure 
prevents  los3  of  cloud  mass  from  be  grid  domain.  However,  clouds  are  observed 
to  be  advected  in  a  direction  different  than  the  streamflow. 


Trajectories  at  the  lower  boundary  are  also  modified  to  account  for 
terrain  effects  on  cloud  formation.  As  described  in  Section  3.4,  the  gradient 
level  trajectories  are  constrained  to  follow  the  terrain  surface.  Extremely 
large  values  of  the  vertical  trajectory  component,  ^p,  are  induced  by  steep 
terrain  features.  This  effect  causes  copious  amounts  of  precipitation  in 
regions  of  steep  ascent  and  very  dry  conditions  in  regions  of  steep  descent. 
This  dry  bia3  on  the  lee  side  of  mountains  is  advected  to  other  regions. 
However,  it  cannot  be  modified  by  a  moisture  source  region  since  no  mechanism 
of  this  type  exists  In  5LAY8R.  To  reduce  these  terrain  induced  effects, 
horizontal  trajectories  are  halved  if  the  vertical  component  exceeds  50  mb/3 
hr.  A  new,  reduced  vertical  component  is  computed  from  these  shortened 
trajectories . 


Forecasters  should  be  aware  that  cloud  formation  and  advection  near  the 
terrain  surface  involves  complex  physical  relationships.  Parameterizing  these 
relationships  produces  less  than  perfect  cloud  forecasts.  However,  in  the 
future,  advanced  computers  will  explicitly  solve  these  relationships, 
resulting  in  Improved  low-level  cloud  forecasts. 


Cloud  and  temperature  advection  is  not  the  only  forecast  mechanism  that 
effects  these  parameters.  Other  sub-scale  mechanisms,  such  as  those  causing 
cumulus  formation  or  the  surface  heating  and  cooling  of  temperatures  over  land 
masses,  are  present  in  the  atmosphere.  Section  4  describes  the  parameteriza¬ 
tion  of  these  and  other  mechanisms.  It  should  be  recognized  that  in  areas 
where  these  effects  are  especially  strong,  cloud  and  temperature  forecasts 
from  5 LAYER  are  not  explicitly  solved.  Once  additional  computer  resources  are 
obtained,  sub-scale  effects  on  clouds  and  temperature  may  be  more  precisely 
defined. 
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Pig.  17.  Depiction  of  how  compaction  of  eighth-mesh  data  to 
half  mesh  data  eliminates  cloud  boundary  distinction.  Eighth-  mesh 
cloud  values  In  percent  are  below  the  grid  points.  Values  above 
the  grid  points  are  half-mesh  cloud  amounts  derived  from 
eighth  mesh  values. 


8 .  SUMMARY 


The  cloud  forecasting  models  now  in  use  at  AFGWC  have  evolved  over  the 
past  20  years  to  meet  Army,  Air  Force,  and  other  Department  of  Defense 
requirements.  Because  clouds  obstruct  visibility,  military  decision  makers 
will  continue  to  need  cloud  forecasts  Increased  under  si and iu&  ef  lit- 
physical  processes  which  cause  clouds  to  form  and  dissipate  as  wt-ll  as 
increased  computer  capacity,  will  provide  greater  eloud  forecasting  skill  In 
the  future. 

The  first  step  in  the  cloud  forecasting  model  is  the  initialization  of  the 
moisture  and  temperature  values  at  the  grid  points.  Moisture  is  derived  from 
the  RTNEPH  and  temperature  is  obtained  from  HIRAS.  Forecasts  are  made  by 
advecting  the  temperature  and  moisture  values  with  3-hour  Lagrangian 
trajectories  from  the  AWS  GSM  model.  Before  the  forecasting  can  begin,  all 
the  variables  must  be  defined  at  the  model  grid  points  by  either  compacting  or 
interpolating.  Currently  the  cloud  data  are  compacted  from  25  nm  to  100  nm 
resolution  and  the  temperature  data  are  interpolated  from  200  nm  to  100  nm 
resolution. 

Forecasts  are  made  by  advecting  air  parcels  along  the  calculated 
trajectories.  The  advected  parcel  temperature  and  moisture  are  modified  by 
vertical  displacement,  entrainment,  evaporation  and  diurnal  effects.  The 
forecast  temperature  and  moisture  are  used  to  derive  forecasts  of  dew-point 
depression,  stability,  precipitation,  cloud  type,  precipitation  type,  and 
icing. 

High-resolution  cloud  forecasts  are  obtained  over  limited  areas  for  short 
forecast  periods  by  interpolating  the  100  nm  resolution  3-hour  moisture 
changes  to  25  nm  resolution  initial  conditions.  These  forecasts  demonstrate  a 
small  but  significant  improvement  in  skill. 

In  the  tropical  regions,  forecasts  of  24-hour  persistence  in  3-hour 
increments  have  out-performed  dynamic  models  except  in  areas  of  moving 
tropical  disturbances. 

The  compilation  of  verification  statistics  for  the  Five  Layer  Model  and 
the  High  Resolution  Cloud  Prognosis  models  document  the  current  forecast 
skill.  The  movement  of  synoptic-scale  systems  depends  on  the  accuracy  of  the 
model  from  which  the  trajectories  are  obtained.  The  parameterizations  of 
entrainment,  evaporation,  convective  development,  and  diural  effects  are 
recognized  weaknesses  which  require  the  constant  attention  of  the  user.  The 
formation  and  dissipation  of  low  clouds  represents  another  major  weakness  of 
the  model  because  these  clouds  are  influenced  by  terrain  features  as  well  as 
diurnal  heating  and  cooling.  A  future  rewrite  of  the  5LAYER  model  will  make 
it  global  in  domain,  increase  the  horizontal  resolution,  and  include  better 
parameterizations  to  correct  the  weaknesses  stated  earlier. 
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10.  APPENDIX  A.  CONDENSATION  PRESSURE  SPREAD  DERIVATION 


The  derivation  of  Pg  begins  by  specifying  the  mixing  ratio,  w,  at  any 
pressure  as 


w  -  p.S/P>  H>X 


(A  I 


where 


a(Td  -  273) 

b  +  T,  -  273 
d 


( A2 


This  formulation,  known  as  Teten's  equation.  Is  found  In  Berry,  Bollay,  and 
Beers  (1945).  Values  of  a  and  b  are  known  constants  and  depend  on  whether  the 
process  occurs  over  Ice  or  water.  Td  is  the  dew-point  temperature  In 
Kelvin.  For  the  case  of  saturation  ws  Is  similarly  described  a3 


wg  =  (3.8/r)10Xs 

where 


(A3: 


a(T  -  273) 
s 

Xs  ~  b  +  T  -  273 
s 

T3  13  the  condensation  temperature  In  Kelvin.  Because  of  the  conservative 
nature  of  the  dry  adiabatic  process,  the  values  w  and  wg  of  an  air  parcel 
lifted  to  condensation  are  unchanged.  Therefore 

(3.8/P)  10X  =  (3.8/P  )10Xf5 

s 

reordering  terms  one  can  write: 

log  Pg  =  log  P  +  (x3-x).  (A5) 
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wvV>v/vv>>>>?.v.‘.v; 


Using  Poisson's  equation: 


Ts  =  T  (Pa/P)°-286 

and  substituting  equations  (A2),  (A4)  and  (A6)  into  (A5)  the  final 
relationship  becomes: 


iog  Pg 


P  0.286 

a[T(-~ )  -  273]  a[Td  -  273] 

P  0.286  "  b  +  T  -  273 

b[T(~)  -  273] 


While  (A7)  cannot  be  directly  solved  for  Pa ,  if  given  known  values  of  P,  T, 
Td,  a,  and  b,  one  can  calculate  approximate  values  of  Pa  in  a  recursive 
fashion.  With  this  approximate  value  of  Pa,  the  variables  in  equation  (3) 
are  now  determined,  and  the  CPS  can  be  calculated. 


11.  APPBNDIX  B.  TABLES  RELATING  CONDENSATION  PRESSURE  SPREAD  AND  CLOUD  AMOUNT 


Table  Bl.  Conversion  of  Cloud  Amount  to  Condensation  Pressure  Spread 


The  following  four  tables  are  used  by  5LAYBR  to  convert  from  cloud  amount  in 
percent  to  CPS  values  in  mb.  Bach  table  is  labeled  with  the  pressure  level  at 
which  it  is  valid.  CPS  values  are  Interior  to  the  tables.  Units  of  cloud 
percent  form  the  abscissa,  while  tens  cf  cloud  percent  appear  on  the  ordinate 
of  each  table. 


Table  Bla.  850  mb 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

120.0 

114.0 

111.0 

108.2 

106.0 

103.2 

100.5 

98.0 

95.7 

93.7 

90.6 

10 

90.6 

88.4 

86.0 

83.0 

81.0 

77.7 

75.5 

73.4 

71.0 

69.6 

68.2 

20 

68.2 

66.8 

66.0 

64.0 

62.6. 

61.0 

59.6 

58.0 

56.6 

55.8 

55.0 

30 

55.0 

53.7 

52.5 

51.6 

50.6 

49.0 

48.2 

47.4 

46.6 

45.8 

45.0 

40 

45.0 

44.2 

43.4 

42.6 

41.8 

41.0 

40.2 

39.4 

38.6 

37.8 

37.0 

50 

37.0 

36.2 

35.4 

34.8 

34-4 

34-1 

33.7 

33.4 

33.0 

32.6 

32.3 

60 

32.3 

31.9 

31.6 

31.2 

30.7 

30.3 

29.7 

29.2 

28.8 

28.3 

27.9 

70 

27.9 

27.5 

27.0 

26.6 

26.1 

25.7 

25.2 

24.7 

24.2 

23.7 

23.2 

80 

23.2 

22.5 

21.9 

21.2 

20.4 

19.5 

18.6 

17.7 

16.8 

15.6 

14.5 

90 

14.5 

13.5 

12.4 

11.3 

10.1 

9.0 

7.6 

6.0 

4.2 

2.6 

1.0 

Table 

Bib. 

700  mb 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

109.0 

107.0 

105.0 

102.8 

100.8 

99.2 

96.7 

95.6 

94.5 

93.4 

92.2 

10 

92.2 

91.0 

90.2 

89.4 

87.0 

86.0 

85.0 

83.2 

81.0 

79.0 

77.0 

20 

77.0 

75.2 

73.3 

71.5 

70.0 

68.3 

66.9 

65.3 

63.9 

62.4 

61.0 

30 

61.0 

59.8 

58.7 

57.5 

56.2 

55.1 

54.0 

52.9 

51.8 

50.7 

49.8 

40 

49.8 

48.8 

48.0 

47.1 

46.2 

45.3 

44.5 

43.7 

42.8 

41.9 

41.0 

50 

41.0 

40.7 

40.3 

40.0 

39.6 

39.3 

38.9 

38.6 

38.2 

37.9 

37.6 

60 

37.6 

37.2 

36.8 

36.4 

36.0 

35.6 

35.2 

34.8 

34.3 

33.6 

33.3 

70 

33.3 

32.8 

32.2 

31.8 

31.3 

30.7 

30.2 

29.6 

29.0 

28.4 

27.9 

80 

27.9 

27.3 

26.7 

26.2 

25.6 

25.0 

24.2 

23.4 

22.4 

21.4 

20.0 

90 

20.0 

18.4 

17.3 

16.0 

14.5 

13.2 

11.4 

9.6 

7.2 

4.7 

1.0 
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Table  Blc.  500  mb 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

101.0 

100.2 

99.4 

98.6 

97.8 

97.1 

96.2 

95.4 

94.5 

93.7 

92.9 

10 

92.9 

92.1 

91.4 

90.6 

89.8 

89.0 

88.0 

87.0 

85.0 

83.0 

81.0 

20 

81.0 

79.3 

77.7 

76.0 

74.3 

72.7 

71.0 

70.2 

69.3 

68.4 

67.5 

30 

67.5 

66.7 

65.8 

64.9 

64.0 

63.2 

62.3 

61.5 

60.0 

58.3 

56.4 

40 

56.4 

54.7 

52.0 

51.0 

50.3 

49.6 

48.9 

48.3 

47.6 

46.9 

46.1 

50 

46.1 

45.5 

44.8 

44.1 

43.4 

42.7 

42.0 

41.3 

40.8 

40.4 

40.0 

60 

40.0 

39.6 

39.2 

38.8 

38.4 

38.0 

37.6 

37.2 

36.8 

36.3 

35.9 

70 

35.9 

35.5 

35.0 

34.6 

34.2 

33.7 

33.3 

32.8 

32.1 

31.5 

30.7 

80 

30.7 

29.7 

28.8 

28.0 

27.2 

26.2 

25.1 

24.0 

23.0 

21.9 

20.7 

90 

20.7 

19.6 

18.3 

16.9 

15.3 

13.4 

11.4 

9.4 

7.4 

5.0 

1.0 

Table 

Bid. 

300  mb 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

97.0 

96.2 

95.4 

94.6 

93.8 

93.2 

92.2 

91.4 

90.5 

89.7 

88.9 

10 

88.9 

88.1 

87.4 

86.6 

85.8 

85.0 

84.5 

84.0 

83.5 

83.0 

81.0 

20 

81.0 

79.3 

77.7 

76.0 

74.3 

72.7 

71.0 

70.2 

69.3 

68.4 

67.5 

30 

67.5 

66.7 

65.8 

64.9 

64.0 

63.2 

62.3 

61.5 

60.0 

58.3 

56.4 

40 

56.4 

54.7 

52.0 

51.0 

50.3 

49.6 

48.9 

48.3 

47.6 

46.9 

46.1 

50 

46.1 

45.5 

44.8 

44.1 

43.4 

42.7 

42.0 

41.3 

40.8 

40.4 

40.0 

60 

40.0 

39.6 

39.2 

38.8 

38.4 

38.0 

37.6 

37.2 

36.8 

36.3 

35.9 

70 

35.9 

35.5 

35.0 

34.6 

34.2 

33.7 

33.3 

32.8 

32.1 

31.5 

30.7 

80 

30.7 

29.7 

28.8 

28.0 

27.2 

26.2 

25.1 

24.0 

23-0 

21.9 

20.7 

90 

20.7 

19.6 

18.3 

16.9 

15.3 

13.4 

11.4 

9.4 

7.4 

5.0 

1.0 
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Table  B2.  Conversion  of  Condensation  Pressure  Spread  to  Percent  of  Cloud 


The  following  four  tables  are  used  by  5LAYBR  to  convert  from  CPS  In  mb  to 
cloud  amount  In  percent.  Bach  table  Is  labeled  with  the  pressure  level  at  which 
It  Is  valid.  Cloud  values  are  Interior  to  the  tables.  Units  of  CPS  form  the 
abscissa,  while  tens  of  CPS  appear  on  the  ordinate  of  each  table. 


Table 

B2a.  850  mb 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9  10 

0 

100.0 

100.0 

99.4 

98.7 

98.1 

97.6 

97.0 

96.4 

95.7 

95.0  94.1 

10 

94.1 

93.3 

92.4 

91.5 

90.5 

89.5 

88.7 

87.8 

86.7 

85.5  84.4 

20 

84.4 

83.3 

81.8 

80.3 

78.4 

76.4 

74.3 

72.1 

69.8 

67.5  65.5 

30 

65.5 

63.5 

60.8 

58.0 

55.2 

52.5 

51.2 

50.0 

48.7 

47.5  46.2 

40 

46.2 

45.0 

43.7 

42.5 

41.2 

40.0 

38.7 

37.5 

36.2 

35.0  34.3 

50 

34.3 

33.7 

32.6 

31.5 

30.8 

30.0 

28.7 

27.5 

27.0 

26.5  25.7 

60 

25.7 

25.0 

24.4 

23.7 

23.0 

22.2 

22.0 

20.8 

20.1 

19.5  18.7 

70 

18.7 

18.0 

17.6 

17.2 

16.7 

16.2 

15.8 

15.3 

14.9 

14.5  14.2 

80 

14.2 

14.0 

13.5 

13.0 

12.7 

12.5 

12.0 

11.6 

11.2 

10.7  10.3 

90 

10.3 

9.8 

9.6 

9.4 

8.8 

8.2 

7.9 

7.5 

7.0 

6.6  6.2 

100 

6.2 

5.8 

5.4 

5.1 

4.7 

4.4 

4.0 

3.6 

3.1 

2.7  2.3 

110 

2.3 

2.0 

1.6 

1.2 

1.0 

.8 

.7 

.5 

.4 

.2  .0 

Table 

B2b .  700  mb 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9  10 

0 

100.0 

100.0 

99.9 

99.7 

99.3 

98.9 

98.5 

98.1 

97.7 

97.3  96.8 

10 

96.8 

96.2 

95.7 

95.2 

94.4 

93.6 

93.0 

92.3 

91.4 

90.5  90.0 

20 

90.0 

89.4 

88.4 

87.5 

86.2 

85.0 

83.3 

81.5 

79.8 

78.0  76.3 

30 

76.3 

74.5 

72.5 

70.5 

68.5 

66.5 

64.0 

61.6 

58.7 

55.8  52.9 

40 

52.9 

50.0 

48.9 

47.8 

46.6 

45.3 

44.2 

43.1 

42.0 

40.8  39.8 

50 

39.8 

38.7 

37.8 

36.9 

36.0 

35.1 

34.2 

33.4 

32.6 

31.7  30.8 

60 

30.8 

30.0 

29.3 

28.6 

27.9 

27.2 

26.6 

25.9 

25.2 

24.6  24.0 

70 

24.0 

23.3 

22.7 

22.2 

21.6 

21.1 

20.5 

20.0 

19.5 

19.0  18.5 

80 

18.5 

18.0 

17.6 

17.1 

16.6 

16.0 

15.0 

14.0 

13.8 

13.5  12.2 

90 

12.2 

11.0 

10.2 

9.4 

8.5 

7.5 

6.6 

5.8 

5.4 

5.1  4.5 

100 

4.5 

3.9 

3.4 

2.9 

2.5 

2.0 

1.5 

1.0 

.5 

.0  .0 

110 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0  .0 
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Table  B2c.  500  mb 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

100.0 

100.0 

99.8 

99.5 

99.3 

99.0 

98.6 

98.2 

97.7 

97.2 

96.7 

10 

96.7 

96.2 

95.7 

95.2 

94.7 

94.2 

93.6 

92.9 

92.2 

91.5 

90.6 

20 

90.6 

89.8 

88.9 

88.0 

87.0 

86.1 

85.2 

84.2 

83.0 

81.7 

80.7 

30 

80.7 

79.7 

78.2 

76.7 

74.4 

72.1 

69.8 

67.5 

65.0 

62.5 

60.0 

40 

60.0 

57.5 

56.0 

54.6 

53.2 

51.7 

50.2 

48.8 

47.4 

45.9 

44.5 

50 

44.5 

43.0 

42.0 

41.9 

41.4 

40.8 

40.2 

39.7 

39.2 

38.6 

38.0 

60 

38.0 

37.5 

36.4 

35.2 

34.0 

32.9 

31.8 

30.6 

29.4 

28.3 

27.2 

70 

27.2 

26.0 

25.4 

24.8 

24.2 

23.6 

23.0 

22.4 

21.8 

21.2 

20.6 

80 

20.6 

20.0 

19.5 

19.0 

18.5 

18.0 

17.5 

17.0 

16.0 

15.0 

13.8 

90 

13.8 

12.5 

11.2 

9.8 

8.6 

7.5 

6.3 

5.1 

3.8 

2.5 

1.2 

00 

1.2 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

10 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

Table  B2d.  300 

mb 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

100.0 

100.0 

99.8 

99.5 

99.3 

99.0 

98.6 

98.2 

97.7 

97.2 

96.7 

10 

96.7 

96.2 

95.7 

95.2 

94.7 

94.2 

93.6 

92.9 

92.2 

91.5 

90.6 

20 

90.6 

89.8 

88.9 

88.0 

87.0 

86.1 

85.2 

84.2 

83.0 

81.7 

80.7 

30 

80.7 

79.7 

78.2 

76.7 

74.4 

72.1 

69.8 

67.5 

65.0 

62.5 

60.0 

40 

60.0 

57.5 

56.0 

54.6 

53.2 

51.7 

50.2 

48.8 

47.4 

45.9 

44.5 

50 

44.5 

43-0 

42.0 

41.9 

41.4 

40.8 

40.2 

39.7 

39.2 

38.6 

38.0 

60 

38.0 

37-5 

36.4 

35.2 

34.0 

32.9 

31.8 

30.6 

29.4 

28.3 

27.2 

70 

27.2 

26.0 

25.4 

24.8 

24.2 

23.6 

23.0 

22.4 

21.8 

21.2 

20.6 

80 

20.6 

20.0 

19.5 

19.0 

17.0 

15.0 

13.8 

12.5 

11.2 

9.8 

8.6 

90 

8.6 

7.5 

6.3 

5.1 

3.8 

2.5 

1.2 

.0 

.0 

.0 

.0 

00 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

10 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 
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12.  APPENDIX  C  WHOLE  KBSH  TO  HALF  HB3H  INTERPOLATION 


The  Interpolation  method  developed  to  Initialise  SLATER  half- mesh  fields 
from  whole-mesh  fields  is  one  where  successive  passes  are  made  on  the  half- amah 
octagon.  Bach  pass  uses  results  of  previous  passes  to  Interpolate  different  sets 
of  points.  The  passes  are  as  follows: 


Note:  <vgt>  denotes  the  point  weight 

#  denotes  a  whole-mesh  point 

♦  denotes  a  previously  Interpolated  half-mesh  point 

•  denotes  point  being  Interpolated 

Pass  #0:  Unpacks  the  whole-mesh  data,  sultlplles  each  point 

by  48  (to  prevent  round-off  errors  with  integer 
arithmetic),  and  assigns  them  to  the  appropriate 
half-mesh  point.  This  pass  Is  performed  once  to  prepare 
the  data  for  Interpolation. 

Pass  #1:  Interpolate  the  diagonal  octagon  boundary  points. 

<l/3>  ♦  ♦  <l/3>  <l/3>  ♦  ♦  <l/3> 


♦  <l/3> 


♦  <l/3> 


Actual  weights  Effective  weights 

Pass  »2:  Interpolate  the  octagon  Interior,  using  only  whole  mesh 

equivalent  points. 


<l/4> 

♦ 

♦  <l/4> 

<l/4>  * 

#  <l/4> 

t 

♦ 

<l/4> 

♦ 

#  <l/4> 

<l/4>  # 

#  <l/4> 

Actual  weights  Effective  weights 

Pass  #3:  Interpolate  remaining  octagon  boundary  points  using  results 

from  Pass  #2. 


<l/3>  #  •  #  <l/3> 

<5/l2>  #  t 

♦  <5/12> 

♦  <l/3> 

♦ 

<1/12>  ♦ 

♦<l/12> 

Aotual  weights  Effective  weights 


Pass  #4: 


<l/4>  # 


Pass  #5: 


Pass  #6: 


Interpolate  along  constant  J  (rows),  between  whole-mesh 
equivalent  points,  using  results  from  Passes  #1  and  #2. 


<l/4> 


<1/16>  # 


#  <1/16> 


*  <l/4> 


<6/l6>  #  @  #  <6/16> 


<l/4> 


Actual  weights 


<1/16>  # 


#  <1/16> 


Effective  weights 


Interpolate  along  constant  I  (columns),  between  whole-mesh 
equivalent  points,  using  results  from  Passes  #1  and  #2. 

<l/4> 


<l/4>  + 


+  <l/4> 


<l/4> 


Actual  weights 


<1/16>  <6/l6> 

*  * 


#  # 

<1/16>  <6/16> 

Effective  weights 


<1/16> 

# 


* 

<1/16> 


Divide  each  half-mesh  point  by  48,  and  pack  into  the 
appropriate  half-mesh  data  array. 


13.  APPENDIX  D.  COMPUTATION  OF  TRAJECTORIES. 


( 


A  derivation  of  the  equations  that  are  used  to  compute  the  trajectory 
displacements  (Section  3)  In  the  5LAYRR  cloud  model  follows. 

Given  an  air  parcel  which  follows  a  three-dimensional  trajectory  In  lt3 
motion  from  a  starting  point  S  to  end  point  E,  we  want  to  express  the  upstream 
displacement  from  point  B  to  point  S  In  an  x,  y,  p  coordinate  system  In  terms 
of  the  mean  wind  components  u,  v,  and  UJ  valid  over  the  time  Interval,  £t,  of 
the  parcel’s  motion:  _ 

Ax  =  -uAt 

Ay  =  -vAt 

(Dl) 

Ap  =  -u)At 

Let  ue,  ve,  and  (Je  be  equal  to  wind  components  at  the  trajectory  end 
point  averaged  over  the  time  Interval,  At.  and  u3,  v3,  and  be  equal 
to  wind  components  at  the  trajectory  starting  point  averaged  over  the  time 
Interval  /it: 


The  wind  components  u3,  v3,  and  t/Jj.  can  be  expressed  In  terms  of  a  Taylor 
series  expansion  In  the  three  coordinate  directions  about  the  point  E  where 
the  terms  containing  partial  derivatives  higher  than  first  order  have  been 
truncated: 


,  3u  .  ,  9u  .  ,  9u  ._ 
u  *  u  +  -r—  Ax  +  x—  Ay  +  AP 
s  e  9x  9y  dP 


v  =  v  +-|^-Ax+~Ay+-i~jAP 


,  9u)  .  .  9(i)  A  , 

(i)  =  (i)  +  -z—  Ax  +  —  Ay  +  AP 

s  e  9x  9y  9P 


The  partial  derivatives  In  equations  (D2)  are  of  the  wind  components 
averaged  over  the  time  Interval  At. 


By  assuming 


u  =  1/2  (u  +  u  ) 
e  s 


v  =  l/2(v  +  u  ) 

e  s 


<i)  =  l/2(u)  +  a)  ) 

e  s 


'  vs'A 


vv 


and  substituting  the  expressions  In  (D2)  for  ua,  va,  and  ,  we  get 

—  1  9u  .  1  9u  .  1  9u 

U  =  ue  +  2  ^  Alt  +  2  37  +  2  3?  AP 

—  ,  1  9v  .  1  9v  A  .  1  9v  . 

v  ‘  ve  +  2  81  Ax  +  2  3?  Ay  +  2  8P  AP  “ 

—  ,  1  9u>  .  .  1  9w  .  ,1  9w 

U)  =  U,e  +  29^Ax  +  29^Ay  +  29PAP 

Equations  (DA)  can  be  written  as  follows  by  substituting  for 
Ax,  Ay.  and  £P  from  equations  (Dl): 

—  1  9u  —  1  9u  .  —  1  9u  — 

u  "  ue  -  2  3i  Atu  '  2  87  Atv  -  2  8?  At" 

—  1  9v  .  -  1  9v  .  —  1  9v  .  —  / r 

V  ■  \  -  2  FF  AtU  -  2  8?  ACV  -  2  8?  At“  (t 

—  1  9u)  .  —  1  9u>  .  —  1  9u>  .  — 

"  ’  “e  ‘  2  aF  Atu  '  2  8^  Atv  -  2  8P  At“ 

Equations  (D5)  form  a  system  of  three  equations  In  three  unknowns,  u,  v,  and 
<J,  which  can  be  written  In  matrix  form  as 


1  9w  , 

1  *iSAt 


1  9V  A 

li?1 2' 


if!  it 

2  9y 


1  A 

2  9P 


Ip-  At 
2  9P 


1  9u)  . 

2  8P  At 


The  matrix  equation  AX  =  B,  has  a  nontrivial  solution  vector  X,  provided  |a|  /  0 
To  get  (D6)  In  a  simplified  form,  we  let 


,  ,  1  9u  .  1  3U 

1  +  2  At  ay  -  2  8j  At  a, 


1  9u 

2  9P  At 


i  4.  I  av  1  9v  .  , 

u!5‘'  ey  ■  2  3?  At  bp 


,  1  3(d 

J  +  -r  —  At 
2  9x 


1  3u  . 

2  8y  At 


1  8u  , 

2  7p  At 


Expanding  the  expressions  for  Dl,  D2,  D3 ,  and  A  we  get 


|  ^ll  ~  ue(eyep  ~  cybp) 
I  a2|  ;  ex<veep  ^ebp> 


-  ay(veep  -  webp)  ap(veCy  -  weey) 
**■  ue(bxep  bpcx)  ♦  ap(bx^e  —  cxV£) 


(Dll) 


A3 1  *x(®t 


CyVe)  ~  ®y (bx  ^  —  cxv^)  +  ue(bxcy  ”  *-x®y ^ 


A  |  ex(eyep  ^-ybp)  -  »y(bxeD  -  cxbD)  ♦  aD(bxCy  -  cxey) 


y '  uxvp  ^xup'  “p'^x'-y  >“xcy  > 


Each  of  the  terms  in  parentheses  in  equations  (Dll)  appears  more  than  once 
hence,  to  simplify  equations  (Dll),  let 


ci  -  ( ®y  ®p  -  Cy  bp) ,  C2  5  (bx  ep  -  cx  bp), 

c3  5  (bx  cy  ~  cx  ®y)»  dl  s  (ve  «p  ”  we  bp> 

d2  *  (ve  Cy  -  <*»e  ey),  d3  s  (bx  ue  -cx  ve). 

Equations  (Dll)  then  become 

I  Ail*-  ue  ct  -  ay  dL  4  ap  d2, 

I A2 | -  ®x  dl  -  ue  c2  ♦  ap  d3, 

I  A3 1 -  -ex  d2  -  ay  dj  ♦  ue  c3, 

|a|  s  ex  c3  -  ay  c2  4  ap  c3. 


Substituting  from  equations  (D13)  into  equations  (D9),  we  get 


_  ueci  -  aydi  apd2 

u  = - 

excl  -  ayc2  +  apc3 

exdl  -  uec2  +  apd3 

v  = - -  ( D1 

exc 1  -  ayc2  *  apc3 

_  - exd2  ~  ayd3  ♦  uec3 

U/= - 

excl  -  ayc2  +  apc3 

Equations  (D7),  (D12),  and  (D14)  together  yield  explicit  solutions  for  u, 
v,  and  c J  which  are  then  used  in  equations  (Dl)  to  compute  the  upstream 
displacements  from  point  e,  a  grid  point,  to  point  8. 

Using  the  forecast  wind  computed  by  a  coarse-mesh  wind  model,  trajectories 
are  computed  explicitly  at  all  coarse-mesh  grid  points  on  the  AFGWC  octagon 
excluding  the  points  on  the  boundaries.  The  goal  is  to  Interpolate  trajectory 
components  from  the  coarse  mesh  to  the  intervening  half-mesh  grid  points. 
However,  before  doing  this,  the  trajectories  computed  explicitly  in  the 
vicinity  of  the  boundaries  must  be  checked  and  modified,  if  necessary,  to 
ensure  that  they  do  not  originate  outside  the  octagon.  (The  coarse-mesh 
trajectories  on  the  boundaries  are  derived  using  the  trajectories  at  grid 
points  indented  one  coarse-mesh  grid  unit  once  the  check  is  accomplished.) 

A  zone  indented  two  coarse-mesh  grid  units  from  the  octagon  boundary  is 
chosen  in  which  the  trajectories  are  to  be  constrained  to  permit  flow  only 
toward  or  parallel  to  the  boundaries.  This  is  easily  done  for  the  rows  and 
columns  of  grid  points  parallel  to  the  x  and  y  coordinate  axes.  The  grid 
points  parallel  to  the  diagonal  boundaries  are  constrained  using  vector 
operations  which  are  described  below.  The  local  unit  vector  perpendicular  to 
each  of  the  diagonal  boundaries  is  defined  as 


u  —  i  0.7071  i  ±  0.7071  j 


(D15 


I 

where  i  and  J  are  unit  vectors  parallel  to  the  x  and  y  coordinate  axes 
respectively.  The  horizontal  component  of  a  three  dimensional  trajectory 
defined  by  the  vector 


D  =  Axl  +  Ayj 

is  examined  to  see  if  it  has  a  component  along  u  directed  Inward. 

Let  0  be  the  angle  between  u  and  D.  By  definition, 

cos  0  - 

1=1  |d| 

-i,  -A 

Consequently,  if  the  vector  dot  product  of  u  and  D  Is  not  zero,  then 
cos  6  Is  not  zero,  and  T)  has  a  component  along  "u  directed  toward 
the  Interior  of  the  grid.  This  component  Is  defined  by  the  vector 


(D16 


(D17 


u  •  D  -*■ 
- -  u 


u  •  u 


The  component  of  D  parallel  to  the  diagonal  boundary  Is  defined  by 


(D18 


M  ■=  D  -  P  (D19) 

The  horizontal  components  of  the  three-dimensional  trajectory  are  then  set 
equal  to  those  of  M.  If  cos  0=0,  then  the  trajectory  is  not  adjusted. 

With  the  above  adjustments  accomplished,  the  trajectories  at  coarse-mesh 
grid  points  on  the  octagon  boundaries  parallel  to  the  x  and  y  axes  are  set 
equal  to  the  trajectories  at  grid  points  Indented  one  coarse-mesh  row. 
Trajectories  at  coarse-mesh  grid  points  on  the  diagonal  boundaries  are 
computed  by  averaging  the  trajectory  components  at  the  two  closest  Indented 
coarse-mesh  grid  points. 

After  the  preceedlng  method  Is  accomplished,  trajectories  are  available  at 
all  coarse-mesh  grid  points.  These  trajectories  are  then  Interpolated  to  all 
of  the  Intervening  half-mesh  grid  points  of  the  APGWC  octagon. 
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